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Abstract. In this survey, we present a geometric description of Lagrangian 

1/^ , and Hamiltonian Mechanics on Lie algebroids. The flexibility of the Lie al- 

C " 3 ' gebroid formalism allows us to analyze systems subject to nonholonomic con- 

C " ' straints, mechanical control systems, Discrete Mechanics and extensions to 

^N) Classical Field Theory within a single framework. Various examples along the 

^...^ , discussion illustrate the soundness of the approach. 
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1. Introduction 

The theory of Lie algebroids and Lie groupoids has proved to be very useful in 
different areas of mathematics including algebraic and differential geometry, alge- 
braic topology, and symmetry analysis. In this survey, we illustrate the wide range 
of applications of this formalism to Mechanics. Specifically, we show how the flexi- 
bility provided by Lie algebroids and groupoids allows us to analyze, within a single 
framework, different classes of situations such as systems subject to nonholonomic 
constraints, mechanical control systems. Discrete Mechanics and Field Theory. 

The notions of Lie algebroid and Lie groupoid allow to study general Lagrangian 
and Hamiltonian systems beyond the ones defined on the tangent and cotangent 
bundles of the configuration manifold, respectively. These include systems deter- 
mined by Lagrangian and Hamiltonian functions defined on Lie algebras. Lie groups, 
Cartesian products of manifolds, and reduced spaces. 

The inclusive feature of the Lie algebroid formalism is particularly relevant for 
the class of Lagrangian systems invariant under the action of a Lie group of symme- 
tries. Given a standard Lagrangian system, one associates to the Lagrangian func- 
tion a Poincare-Cartan symplectic form and an energy function using the particular 
geometry of the tangent bundle. The dynamics is then obtained as the Hamiltonian 
vector field associated to the energy function through the Poincare-Cartan form. 
The reduction by the Lie group action of the dynamics of this system yields a 
reduced dynamics evolving on a quotient space (which is not a tangent bundle). 
However, the interplay between the geometry of this quotient space and the reduced 
dynamics is not as transparent as in the tangent bundle case. Recent efforts have 
lead to a unifying geometric framework to overcome this drawback. It is precisely 
the underlying structure of Lie algebroid on the phase space what allows a unified 
treatment. This idea was introduced by Weinstein (see also EU), who devel- 
oped a generalized theory of Lagrangian Mechanics on Lie algebroids. He obtained 
the equations of motion using the linear Poisson structure on the dual of the Lie 
algebroid and the Legendre transformation associated with a (regular) Lagrangian. 
In [^, Weinstein also posed the question of whether it was possible to develop a 
treatment on Lie algebroids and groupoids similar to Klein's formalism for ordinary 
Lagrangian Mechanics This question was answered positively by E. Martinez 
in PHI (see also ^| 123 ^2 EH ) ■ The main notion was that of prolongation of a Lie 
algebroid over a mapping, introduced by P.J. Higgins and K. Mackenzie ^Hl- More 
recently, the work has developed a description of Hamiltonian and Lagrangian 
dynamics on a Lie algebroid in terms of Lagrangian submanifolds of symplectic Lie 
algebroids. An alternative approach, using the linear Poisson structure on the dual 
of the Lie algebroid was discussed in P^ . 
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In this paper, we present an overview of these developments. We make special 
emphasis on one of the main advantages of the Lie algebroid formalism: the pos- 
sibility of establishing appropriate maps (called morphisms) between systems that 
respect the structure of the phase space, and allow to relate their respective prop- 
erties. As we will show, this will allow us to present a comprehensive study of the 
reduction process of Lagrangian systems while remaining within the same category 
of mathematical objects. 

Wc also consider nonholonomic systems (i.e., systems subject to constraints in- 
volving the velocities, see ^U] for a list of references) and control systems evolving 
on Lie algebroids ^Hl- This is motivated by the renewed interest in the study of 
nonholonomic mechanical systems for new applications in the areas of robotics and 
control. In particular, we provide widely applicable tests to decide the accessi- 
bility and controllability properties of mechanical control systems defined on Lie 
algebroids. 

We end this survey by paying attention to two recent developments in the context 
of Lie algebroids and groupoids: Discrete Mechanics and Classical Field Theory. 
Discrete Mechanics seeks to develop a complete discrete-time counterpart of the 
usual continuous-time treatment of Mechanics. The ultimate objective of this effort 
is the construction of numerical integrators for Lagrangian and Hamiltonian systems 
(see |37] and references therein). Up to now, this effort has been mainly focused 
on the case of discrete Lagrangian functions defined on the Cartesian product of 
the configuration manifold with itself. This Cartesian product is just an example 
of a Lie groupoid. Here, wc review the recent developments in |85| . where we 
proposed a complete description of Lagrangian and Hamiltonian Mechanics on Lie 
groupoids. In particular, this description covers the analysis of discrete systems 
with symmetries, and naturally produces reduced geometric integrators. Another 
extension that we consider is the study of Classical Field Theory on Lie algebroids. 
Thinking of a Lie algebroid as a substitute of the tangent bundle of a manifold, we 
substitute the classical notion of fibration bundle by a surjective morphism of Lie 
algebroids tt : E — > F. Then, we construct the jet space Jtt as the affine bundle 
whose elements are linear maps from a fiber of to a fiber of E, i.e., sections of 
the projection tt. After a suitable choice of the space of variations, we derive the 
Eulcr-Lagrange equations for this problem. 

The paper is organized as follows. In Section |21 we present some basic facts on 
Lie algebroids, including results from differential calculus, morphisms and prolon- 
gations of Lie algebroids, and linear connections. Wc also introduce the notion 
of Lie groupoid; Lie algebroid associated to a Lie groupoid, and morphisms and 
prolongations of Lie groupoids. Various examples arc given in order to illustrate 
the generality of the theory. In Section O we give a brief introduction to the La- 
grangian formalism of Mechanics on Lie algebroids, determined by a Lagrangian 
function L : E — > R on the Lie algebroid r : E — > M. Likewise, we introduce the 
Hamiltonian formalism on Lie algebroids, determined by a Hamiltonian function 
H: E* — > R, where t* : E* — > M is the dual of the Lie algebroid E — > M. In 
Section ^ we introduce the class of nonholonomic Lagrangian systems. We study 
the existence and uniqueness of solutions, and characterize the notion of regularity 
of a nonholonomic system. Under this property, wc derive a procedure to obtain 
the solution of the nonholonomic problem from the solution of the free problem by 
means of projection techniques. Moreover, we construct a nonholonomic bracket 
that measures the evolution of the obscrvablcs. and wc study the reduction of 
nonholonomic systems in terms of morphisms of Lie algebroids. In Section [31 we 
introduce the class of mechanical control systems defined on a Lie algebroid. We 



4 J. CORTES, M. DE LEON, J. C. MARRERO, D. MARTfN DE DIEGO, AND E. MARTfNEZ 



generalize the notion of affinc connection control system to the setting of Lie alge- 
broids, and introduce the notions of (base) accessibility and (base) controllability. 
We provide sufficient conditions to check these properties for a given mechanical 
control system. In Section we study Discrete Mechanics on Lie groupoids. In 
particular, we construct the discrete Euler-Lagrange equations, discrete Poincare- 
Cartan sections, discrete Lcgendrc transformations, and Noethcr's theorem, and 
identify the preservation properties of the discrete flow. In the last section, we 
extend the variational formalism for Classical Field Theory to the setting of Lie 
algebroids. Given a Lagrangian function, we study the problem of finding criti- 
cal points of the action functional when we restrict the fields to be morphisms of 
Lie algebroids. Throughout the paper, various examples illustrate the results. We 
conclude the paper by identifying future directions of research. 



2. Lie algebroids and Lie groupoids 



2.1. Lie algebroids. Given a real vector bundle t: E — > AI, let Scc(t) denote 
the space of the global cross sections oi t: E — > M . A Lie algebroid E over a 
manifold M is a real vector bundle r: E — > M together with a Lie bracket |-, ■] 
on Sec(T) and a bundle map p: E — > TM over the identity, called the anchor 
map, such that the homomorphism (denoted also p: Sec(T) — > X(M)) of C°°(A/)- 
modules induced by the anchor map verifies 

lXJY\^flX,Yl+p{X){J)Y, 

for X,Y & Sec(T) and / £ C°°(M). The triple {E, |-, •], p) is called a Lie algebroid 
over M (see |21||211). If {E, |-, •], p) is a Lie algebroid over M, then the anchor map 
p: Sec(T) — > X{M) is a homomorphism between the Lie algebras (Sec(r), |-, •]) 
and (X(M), [•,•]). 

In what concerns to Mechanics, it is convenient to think of a Lie algebroid as 
a generalization of the tangent bundle of M . One regards an clement a of as 
a generalized velocity, and the actual velocity v is obtained when applying the 
anchor to a, i.e., v = p{a). A curve a: [toj^i] — * E is said to be admissible if 
mit) = p{a{t)), where ■m{t) = r(a(t)) is the base curve. 

Given local coordinates (a;*) in the base manifold M and a local basis of sections 
(bq,) of E, then local coordinates of a point a € E are (x*, y") where a = y°'ea{T{a)). 
In local form, the Lie algebroid structure is determined by the local functions 
and C^^ on M. Both are determined by the relations 

, d 



P(ea) ^ 

and they satisfy the following equations 



Pa 



dp] 



rfi P _ r? ' — r 



dx^ 



and 



E 



cyclic(a,/3,7) 



Pa 



= 0. 



(2.1) 
(2.2) 

(2.3) 



Cartan calculus. One may define the exterior differential of E, d: Sec(A'^T*) 
Scc(a''+^t*), as follows 



±oiXo, . . . , Xfe) = ^(-l)V(^O(^(^o, . . . , ^z, . . . , Xk)) 



i=0 



Y,i-^Y^=^ilX,,X,lXo,...,Xi,---Xj,---,Xk), 



(2.4) 
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for iu G Sec(A'"'r*) and Xq, . . . ,Xk £ Sec(T). c? is a cohomology operator, that is, 
(P = 0. In particular, if / : M — > M is a real smooth function then df{X) = p{X)f, 
for X e Sec(T). Locally, 

dx' = p^e" and de"^ = -^C-^/je" A e'^, 

where {e"} is the dual basis of {ea}- We may also define the Lie derivative with 
respect to a section X of i? as the operator Lx- Sec(A'^T*) — > Sec(A'^r*) given 
by l^x = ix ° d + d o ix (for more details, see |32l 

Morphisms. Let {E, |-, (resp., {E' , |-, ■]' ,p'j) be a Lie algebroid over a man- 
ifold AI (resp., M') and suppose that E — > E' is a vector bundle morphism 
over the map 'fo ■ — > A/'. Then, the pair (^P, ^fo) is said to be a Lie algebroid 
morphism if 

fi((1', *o)*0') = (*, *o)*(rf'0'), for all cj)' € Sec(A'= (£')*) and for all fc, (2.5) 

where d (resp., d') is the differential of the Lie algebroid E (resp., E') (see UHl)- 
Note that {^,^o)*(p' is the section of the vector bundle a''E* — > M defined for 
> by 

{{^, 'fo)*(f>')x{ai, . . . , flfc) = <?!''3,„(^)(*(ai), . . . , ^'(afc)), 

for X e M and ai, . . . , ak e E^, and by (*, '^o)*f = / o *o for / G Scc(a"£"*) = 
C°°(A/'). In the particular case when M = M' and *o = then (|2.5|) holds if 
and only if 

|*oX,*or]' = 5'|X,r], p'{^X)^ p{X), for X,y G Scc(t). 

Linear connections on Lie algebroids. Let r : E — > A/ be a Lie algebroid over 
A/. A connection on E \s a, M-bilinear map V: Sec(i?) x Sec(£') — > Sec(i?) such 
that 

VfxY^fVxY, VxifY)^ p{X){f)Y + fVxY 
for / G C°°(Af ) and X, F G Scc(£;). 

Given a local basis {cq} of Scc(i?) such that X = X^e^ and F = y'e^ then 

^xY = x- (^p^^ + rl,y/^^ e,. 

The terms F^^ are called the connection coefficients. The symmetric product 
associated with V is given by 

{X -.Y) =yxY + VyX, X,YeSec{E). 

Since the connection is C°° (Af )-linear in the first argument, it is possible to 
define the derivative of a section Y G Sec(i?) with respect to an element a G E,n 
by simply putting 

VaY = {VxY)im), 

with X G Sec(i?) satisfying X(m) — a. Moreover, the connection allows us to 
take the derivative of sections along maps and, as a particular case, of sections 
along curves. If we have a morphism of Lie algebroids $ : F — > E over the map 
ip : N — > M and a section X : N — > E along ip (i.e., X{ti) G for n G N), 

then X may be written as 

p 

X-5^F,(X;0^), 

1 = 1 
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for some sections {^i, . . . , Xp} of E and for some functions Fi, . . . ,Fp G C°°{N), 
and the derivative of X along (p is given by 
p 

VbX = [iPF{b)Fi)Xiiipin)) + Fi(n)V$(b)Xi] , for & G F„ , 

1=1 

where pp is the anchor map of the Lie algebroid F — > N. 

A particular case of the above general situation is the following. Let a : I — > E 
be an admissible curve and b : / — > E he a curve in E, both of them projecting by 
T onto the same base curve in M, T{a{t)) = m{t) ~ T{h{t)). Take the Lie algebroid 
structure TI — > I and consider the morphism <I> : TI — > E, ^{t,i) = ia{t) over 
m : I — > M. Then one can define the derivative of b{t) along a{t) as '^d/dtb{t)- 
This derivative is usually denoted by Va{t)b{t). In local coordinates, this reads 

e-y(m(t)), for all t. 

The admissible curve a : I — > E is said to be a geodesic for V if Va{t)a{t) = 
(see [13). 

Now, let 5 '■ E Xm E — > M be a bundle metric on a Lie algebroid r : E — > M. 
In a parallel way to the situation in the tangent bundle geometry, one can see that 
there is a canonical connection on E associated with S- In fact, the connection 
is determined by the formula 

2g(V|r,Z) = p{X){5iY,Z))+p{Y)i9{X,Z))-piZ){SiX,Y)) 

+9ix, iz, ¥]) + g(y, |z, xj) - g(z, |f, xj), 

for X,Y,Ze Sec{E). is a torsion-less connection and it is metric with respect 
to g. In other words 

lx,Yj = vsy-v^x, 
pix){9{Y,z)) = g(v|r,z) + g(y,vsz). 

is called the Levi-Civita connection of g (see 

Finally, suppose that E = D (B D'^, with D and D'^ vector subbundles of E, 
and denote by P : i? — > D and Q : E — > D'^ the corresponding complementary 
projectors induced by the decomposition. Then, the constrained connection is 
the connection V on i? defined by 

WxY^P(VxY) + Vx{QY), 

for X,Y G Sec{E) (for the properties of the constrained connection V, see 
Examples. We will present some examples of Lie algebroids. 

1. - Real Lie algebras of finite dimension. Let g be a real Lie algebra of 
finite dimension. Then, it is clear that g is a Lie algebroid over a single point. 

2. - The tangent bundle. Let TM be the tangent bundle of a manifold AI. 
Then, the triple {TM, [■,■], idxM) is a Lie algebroid over A/, where idxAi ■ TM — > 
TM is the identity map. 

3. - Foliations. Let 5" be a foliation of finite dimension on a manifold P and 
Tg^ : — > P be the tangent bundle to the foliation 7. Then, rgr : TJ — > P is a 
Lie algebroid over P. The anchor map is the canonical inclusion : — > TP 
and the Lie bracket on the space Sec(Tg:^) is the restriction to Sec(Tgr) of the standard 
Lie bracket of vector fields on P. In particular, if tt : P — > M is a fibration, 
Tp : TP — > P is the canonical projection and (Tp)|y^ : Vn — > P is the restriction 
of Tp to the vertical bundle to tt, then (Tp)[y^ : Vtt — > P is a Lie algebroid over P. 



^a(t)b{t) 



dt 
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4. - Atiyah algebroids. Let p : Q — > M be a principal G-bundle. Denote by 
$ : G X Q — > Q the free action of G on Q and by T$ : G x TQ — > TQ the 
tangent action of G on TQ. Then, one may consider the quotient vector bundle 
tq\G : TQ/G — > M = Q/G, and the sections of this vector bundle may be 
identified with the vector fields on Q which are invariant under the action <f>. Using 
that every G-invariant vector field on Q is p-projectablc and that the usual Lie 
bracket on vector fields is closed with respect to G-invariant vector fields, we can 
induce a Lie algebroid structure on TQ/G. This Lie algebroid is called the Atiyah 
algebroid associated with the principal G-bundle p : Q — > M (see |26l IH2j 'l. 

5. - Action Lie algebroids. Let {E, |-, -J, p) be a Lie algebroid over a manifold 
M and /: M' — > M be a smooth map. Then, the pull-back of E over /, f*E = 
{ (a;', a) e M' x E \ f{x') = T(a) } , is a vector bundle over M' whose vector bundle 
projection is the restriction to f*E of the first canonical projection pri : M' x E ^ 
M' . However, f*E is not, in general, a Lie algebroid. 

Now, suppose that Sec(T) — > X(M') is an action of E on /, that is, $ is a 
R-linear map which satisfies the following conditions 

^{hx)^{ho f)<^x, = [$x,$r], ^x{hof)^p{x){h)of, 

for X, y G Sec(T) and h G G°°(M). Then, one may introduce a Lie algebroid 
structure (|-, •]$, p$) on the vector bundle f*E M' which is characterized by the 
following conditions 

|Xo/,yo/]<j, = |X,r]o/, p4,(Xo/) = $(X), for X,rGSec(T). (2.6) 

The resultant Lie algebroid is denoted by E t< f and we call it an action Lie 
algebroid (for more details, see |2fijV 

6. - The prolongation of a Lie algebroid over a fibration ^USHlEni- Let 

{E, |-, •], p) be a Lie algebroid over a manifold M and tt : P — > M be a fibration. 
We consider the subset of E x TP 

TfP = { (6, G X TpP I p{b) = Tp7r{v) } , 

where Ttt : TP — > TM is the tangent map to tt, p G Px and 7r(p) = x . We will 
frequently use the redundant notation {p, b, v) to denote the element (6, v) G T^P. 
T^P = Up^zpTjfP is a vector bundle over P and the vector bundle projection Tp 
is just the projection onto the first factor. The anchor of T^P is the projection 
onto the third factor, that is, the map p'^ : T^P — > TP given by p^{p, b, v) = v. 
The projection onto the second factor will be denoted by Ttt : T^P — > E, and it 
is a morphism of Lie algebroids over tt. Explicitly, Ttt{p, b, v) = b. 

An element z G T^P is said to be vertical if it projects to zero, that is Ttt{z) = 
0. Therefore it is of the form [p, 0, v), with v a tt- vertical vector tangent to P at p. 

Given local coordinates (a;*, u^) on P and a local basis {cq} of sections of E, we 
can define a local basis {Xq, V^} of sections of T^P by 



3Ca(p) = (p,ea(7r(p)),pj,^ ^ and Va(p) = 0, 



d 



li z — {p,b,v) is an element of T^P, with b = z^Eq., then v is of the form v 
PlvZ^-r— + V^jAt-, and we can write 

z = z^X^{p) + v^Va{p). 



Vertical elements are linear combinations of {V^}. 
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The anchor map apphed to a section Z of P with local expression Z = 
Z"Xq, + V^Va is the vector field on P whose coordinate expression is 

^^(^)=^^"^"& + ^'^- 

Next, we will see that it is possible to induce a Lie bracket structure on the 
space of sections of T^P. For that, we say that a section X oi Tp : T^P — *■ P 
is projectable if there exists a section X t : E — > AI and a vector field 
U € X{P) which is 7r-projectable to the vector field p{X) and such that X{p) = 
{X{tt{p)), U{p)), for all p G P. For such a projectable section X, we will use the 
following notation X = {X, U) . It is easy to prove that one may choose a local 
basis of projectable sections of the space Sec(Tp). 

The Lie bracket of two projectable sections Zi ~ (Xi, Ui) and Z2 = (X2, U2) is 
then given by 

lZi,Z2Y{p)^[pAXi,X2l{x),[UuU2][p)). P€P, x = 7r{p). 

Since any section of T^P can be locally written as a linear combination of pro- 
jectable sections, the definition of the Lie bracket for arbitrary sections of T^P 
follows. In particular, the Lie brackets of the elements of the basis are 

IX„,X;3r = C2^X^, IX^,VbV^O and lK4,VBf-0, 
and, therefore, the exterior differential is determined by 

dX'^ = -ic^^X" AX/^, dV^ = 0, 

where {X", V"^} is the dual basis to {Xq, V^i}. 

The Lie algebroid T^P is called the prolongation of E over tt or the E- 
tangent bundle to n. 

2.2. Lie groupoids. In this section, we review the definition of a Lie groupoid and 
present some basic facts generalities about them (see OHI for more details) . A 
groupoid over a set M is a set G together with the following structural maps: 

• A pair of maps a : G — > M, the source, and (3 : G — > M, the target. 

These maps define the set of composablc pairs 

G2^{{9,h)eGxG \ f3{g)^a{h)}. 

• A multiplication m : G2 — > G, to be denoted simply by m{g, h) — gh, 
such that 

— a{gh) — a{g) and P{gh) ~ f3{h), 

— g{hk) ~ {gh)k. 

• An identity section e : M — > G such that 

- <^{a{g))g ^ g and ge{P{g)) ^ g. 

• An inversion map i : G — > G, to be denoted simply by i{g) = g^^, 
such that 

- .9"\9 = e(/3(.9)) and gg'^ = e{a{g)). 

A groupoid G over a set M will be denoted simply by the symbol G =t M. 

The groupoid G =t M is said to be a Lie groupoid if G and M are manifolds and 
all the structural maps are differentiable with a and /3 differentiable submersions. 
If G =^ iW is a Lie groupoid then m is a submersion, e is an immersion and i is 
a diffeomorphism. Moreover, if x G M, a~'^{x) (resp., P~^{x)) will be said the 
a-fiber (resp., the (3-fiber) of x. 
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On the other hand, if 5 £ G then the left-translation by g <E G and the 
right-translation by g are the difFeomorphisms 

Ig : a-^{/3{g)) — >a-^{a{g)) ; h — > lg{h) ^ gh, 
rg:(3-\a{g))—.p-\p{g)) ; h—.rg{h) = hg. 

Note that l^"^ ~ Ig-i and r^^ = r^-i. 

A vector field X on G is said to be left-invariant (resp., right-invariant) if 

it is tangent to the fibers of a (resp.. (3) and X(gh) = {Thlg){Xh) (resp., X{gh) = 
(V;,) (1(g)), for ig,h)eG2. 

Now, we will recall the definition of the Lie algebroid associated with G. 

We consider the vector bundle t : Eq — > M, whose fiber at a point x £ M 
is [Eg)x = Vt(x)(^ = Ker{T^(^^'^a). It is easy to prove that there exists a bijection 
between the space Sec(T) and the set of left-invariant (resp., right-invariant) vector 
fields on G. If X is a section of t : Eq — > A/, the corresponding left-invariant 
(resp., right-invariant) vector field on G will be denoted X (resp., X), where 

X{g) = {T,(p(g))lg){X{fi{g))), (2.7) 

X{g) = -(T,(„(,))r,)((T,(„(,))z)(X(a(5)))), (2.8) 
for g G G. Using the above facts, we may introduce a Lie algebroid structure 
(|-, •],p) on Eg, which is defined by 

\Xyi^[X,Y], p{X){x) = {T,(,)p){X{x)), (2.9) 
for X,Y £ Sec(r) and x & M (for more details, see [HlESl)- 

Given two Lie groupoids G ^ M and G' ^ M', a morphism of Lie groupoids 
is a smooth map ^' : G — > G' such that 

{g,h)eG2=^{^ig),^{h))e{G')2 

and 

^(gh) = ^{g)^ih). 

A morphism of Lie groupoids 5" : G — > G' induces a smooth map $0 : M — > M' 
in such a way that 

a' o v]/ = $0 o a, /3' o = $0 o /?, vfi o e = e' o c&o, 

a, /3 and e (resp., a', /?' and e') being the source, the target and the identity section 
of G (resp., G'). 

Suppose that (^P, $0) is a morphism between the Lie groupoids G ^ M and 
G' ^ M' and that t : Eq — > M (resp., r' : Ec — > A/') is the Lie algebroid of 
G (resp., G'). Then, if x G A/ we may consider the linear map ■ {Eg)x — > 
{Eg')<s>o(x) defined by 

'^>xive{x)) = (T,(,)«')(i;,(,)), for G A,G. (2.10) 

In fact, we have that the pair ($, <i>o) is a morphism between the Lie algebroids 
t:Eg — > M and r' : Eg' — > M' (see 1221231). 
Examples. We will present some examples of Lie groupoids. 

1. - Lie groups. Any Lie group G is a Lie groupoid over {e}, the identity element 
of G. The Lie algebroid associated with G is just the Lie algebra q of G. 

2. - The pair or banal groupoid. Let Af be a manifold. The product manifold 
M X M is a Lie groupoid over M in the following way: a is the projection onto 
the first factor and (5 is the projection onto the second factor; e{x) = {x, x), for all 
X G M, ■m{{x,y),{y,z)) = {x,z), for {x,y),{y,z) e M x M and i{x,y) = {y,x). 
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M X Af ^ M is called the pair or banal groupoid. If a; is a point of M, it follows 
that 

= {OJ X T,Af C T^M X T^M ^ r(,,,)(Af x M). 
Thus, the linear maps 

■ T^M — > Ve(^)a, Vx — > {Ox,Vj;), 

induce an isomorphism (over the identity of M) between the Lie algebroids tm '■ 
TM — > M and r : Emxm — > M. 

3. - The Lie groupoid associated with a fibration. Let tt : P > M be 
a fibration, that is, tt is a surjective submersion and denote by G-^ the subset of 
P X P given by 

G^ = {{p,p')ePxP \n{p)^n{p')}. 

Then, is a Lie groupoid over P and the structural maps aj^, p^^, m^j, e^r and 
are the restrictions to G^r of the structural maps of the pair groupoid P x P ^ P. 

If p is a point of P it follows that 

K,(p)a. = { (Op, Yp) e TpP X TpP I {Tp7r){Yp) = } . 

Thus, if (Tp)|y^ : Vt: — > P is the vertical bundle to tt then the linear maps 

($^)p : VpTT — > V;^(p)Q;^, Yp — > (Op, Yp) 

induce an isomorphism (over the identity of M) between the Lie algebroids (Tp)|y^ : 
Vtt > P and r : Eq^ > P- 

4. - Atiyah or gauge groupoids. Let p : Q — > M be a principal left G-bundle. 
Then, the free action $ : G x Q — > Q, {g, q) — > $(g, q) = gq, of G on Q induces, 
in a natural way, a free action $ x $ : G x (Q x Q) — > Q x Q of G on Q x Q given 
by ($ X 3>)(g, (g, q')) = {gq, gq'), for g e G and (q, q') e Q x Q. Moreover, one may 
consider the quotient manifold {Q x Q)/G which admits a Lie groupoid structure 
over M with structural maps given by 



d:{Qx Q)/G — > M 

/3 : (Q X Q)/G M 

e : M ^ (Q X Q)/G 

m:{{Qx Q)/G)2 {Q x Q)/G 

i:{Qx Q)/G -^{Qx Q)/G 



[{q,q')]^piq), 
[{q,q')]^p{q'), 
X — >[iq,q)], ifp('?) = a;, 

[{q,q')]^[{q',q)]. 



This Lie groupoid is called the Atiyah (gauge) groupoid associated with the 
principal G-bundle p : Q — > M (see 

If a; is a point of M such that p{q) = x, with q ^ Q, and pqxq : Q x Q — > 
{Q X Q)/G is the canonical projection then it is clear that 

Vg(:r)a = iT{q,q)PQxQ){{Og} X TgQ) . 

Thus, if tq|G : TQ/G — > M is the Atiyah algebroid associated with the principal 
G-bundle p : G — > M then the linear maps 

{TQ/G)x — > Vi(x)d ; [vg] — > {T(^q,q)PQxQ){Oq,Vg), with Vq e TqQ, 

induce an isomorphism (over the identity of M) between the Lie algebroids r : 
EiQxQ)/G ^ M and tq\G : TQ/G M. 

5.- Action Lie groupoids. Let G ^ M be a Lie groupoid and / : M' — > M 
be a smooth map. If M' f x „ G ^ { {p, g) & P x G \ f{p) = a{g) } then a right 
action of G on / is a smooth map 

M' jX^G^M', (x',g)—>x'g. 
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which satisfies the foUowing relations 

fix'g) = Pig), for (.T',g) eM'/X„G, 

(x'g)h = x'{gh), for {x',g) £ M' /X^ G and {g,h) e G2, and 
x'e{f{x')) = x', for x' e M' . 

Given such an action one constructs the action Lie groupoid M' fX^G over 
M' by defining 



5/ : M' fXa G — >M' 

Pf : Af /x„ G — > M' 

if : M' — > M' /x„ G 

rhf : (Af /x„ G)2 — > M' /X„ G 

i/ : M' fXa G — > M' /x^ G 

Now, if x' G A/', we consider the map x' 



{x',g) — ^ x', 

{x',g) — > x'g, 

x' ~^{x',e{f{x'))), 

{{x',g),ix'g,h))^ix\gh), 

{x',g) — > {x'g,g^^). 



a~Hf{x')) > given by 

x' ■ (5) = x'g. 

Then, if r : Eq — > M is the Lie algebroid of G, the K-finear map <f> : Sec(T) — > 
X(Af ') defined by 

<P{X)ix') ^ (re(/(.'))^' -K^ifi^')))^ for ^ e Sec(T) and x' e A/', 

induces an action of Eg on / : A/' — > A/. In addition, the Lie algebroid associated 
with the Lie groupoid M' fX a G ^ P is the action Lie algebroid Eg x f (for more 
details, see [H^^. 

6.- The prolongation of a Lie groupoid over a fibration. Given a Lie 
groupoid G =^ A/ and a fibration tt : P — > Af, we consider the set 

P X P = P ,x„ G 0X, P = { e P X G X P I ^(p) = Pig) = 7r(p') } . 

G . . . 

Then, P x P is a Lie groupoid over P with structural maps given by 

G 



d" -.P X P — > P 
p-^ : Px P — > P 
^TT . p — , px P 

G G 

m'^ : iP X P)2 — > P X P 

G G 

i'^ : P X P — > P X P 



iPidiP') — 'Pi 
iP,9,P') — 'P\ 
P — > (P,e(7r(p)),p), 
iip,g,p'),ip',h,p")) ip,gh,p"), 

ip,g,p') — ' ip',9^^,p)- 



P X P is called the prolongation of G over tt : P 



M. 



Now, denote by t : Eg — > A/ the Lie algebroid of G, by g the Lie algebroid 



PxP 

M over the fibration tt. If 



of P X P and by T^'^ P the prolongation of t : Eg 
p E P and m = 7r(p), then it follows that 

(eo) - { (Op,i;,(™),Xp) e TpP X iEG)m X TpP \ iTp7r)iXp) = (re(,„)/?)(«e(™)) } 

and, thus, one may consider the linear isomorphism 

(<f>^)p: (£;^g^)p^r/°P, (Op,t;,(„),Xp)^ Xp). (2.11) 



In addition, one may prove that the maps i^^)p, p £ P, induce an isomorphism 
: P G — > r^op between the Lie algebroids E g and T^^p ('for more 

PxP PxP 

details, see [T^'l. 
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A particular case. Next, suppose that P = Eq and that the map tt is just the 
vector bundle projection t : Eg — > M. In this case, 

G 

Eg X Eg = Eg r'Xa G pXr Eg 

G 

and we may define the map 8 : Eg x Eg — > VP (Bg Vet as foUows 

0(We(a(g)),5,'fe(/3(g))) = iiTe{a(g)) {r g O «)) (w£(a(g)) ) , (^e(/3(g)) ^s) (^'e(/3(g)) ) ) , 

for (ue(a(g)),5,W£(/3(g))) £ (£^G)a(g) X G X (£;G)/3(g)- 6 is a bijective map and 

G 

for {Xg,Yg) e V"g/3 © Vg . Thus, the spaces Eg x and ©g T^a may be 
identified and, under this identification, the structural maps of the Lie groupoid 
structure on V/3 <S)g Vol are given by 



{Xg.Yg) ^ {Tglg-^){Yg), 

Veix) > ((r'e(j,)i)(«e(2;)),W£(2;)), 

(x„y,)^((r,z)(r,),(T,i)(x,)), 



a" : VP ©G Va — > Eg 

: VP ©G Va — > Eg 

-.Eg — > VP ©G Va 
r : Vp ©G Va — > Vp ©G Va 

and the multiplication m'^ : {VP ©g Va)2 — > ©g Va is 

m-((X„y,),((^g(rg,o^))(y,),Z,0) = ((rgr,)(Xg),(T,Z,)(Z,0). 

This Lie groupoid structure was considered by Saunders |53| . We remark that the 
G 

Lie algebroid of Eg x Eg = ©g V« =| Sg is isomorphic to the prolongation 
T^^Eg of £;g over t : Eg — > M. 



3. Mechanics on Lie algebroids 



We recall that a symplectic section on a vector bundle tt : F — > M is a section 
UJ of A^TT* which is regular at every point when it is considered as a bilinear form. 
By a symplectic Lie algebroid we mean a pair {E, uj) where r : E — > M is a 
Lie algebroid and a; is a symplectic section on the vector bundle E satisfying the 
compatibility condition dio = 0, where d is the exterior differential of E. 

On a symplectic Lie algebroid {E, uj) we can define a dynamical system for every 
function on the base, as in the standard case of a tangent bundle. Given a function 
H S C°°{M) there is a unique section (Jh S Sec(T) such that 

The section (Jh is said to be the Hamiltonian section defined by H and the 
vector field Xh = p(o'h) is said to be the Hamiltonian vector field defined by 
H. In this way we get the dynamical system x = Xh{x). 

A symplectic structure w on a Lie algebroid E defines a Poisson bracket { , 
on the base manifold M as follows. Given two functions F,G E C°°{M) we define 
the bracket 

{F,Gr =Lo{aF,aG). 
It is easy to see that the closure condition du — implies that { , is a Poisson 
structure on M . In other words, if we denote by A the inverse of to as bilinear form, 
then {F,G}'^ = A{dF,dG). The Hamiltonian dynamical system associated to H 
can be written in terms of the Poisson bracket as i = {x, H}^ . 

By a symplectomorphism between two symplectic Lie algebroids {E,lj) and 
{E' ,Lu') we mean an isomorphism of Lie algebroids (4*, ^'o) from E to E' such that 
^t))*uj' = u!. In this case the base map ^'o is a Poisson diffeomorphism, that is, 
it satisfies '^*o{F' ,G'}'^' = {^-JF', ^^G'}", for aU F',G' e C°°{M'). 
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Sections 13.11 and 13.21 describe two particular and important cases of the above 
construction. 

3.1. Lagrangian Mechanics. In (see also a geometric formalism for 

Lagrangian Mechanics on Lie algebroids was introduced. It is developed in the 
prolongation T^E of a Lie algebroid E over the vector bundle projection t ; E — > 
M. The canonical geometrical structures defined on T^E are the following: 

• The vertical lift C ■ t* E — > T'^E given by {a, b) = (a, 0, where 

is the vector tangent to the curve a + <6 at t = 0. 

• The vertical endomorphism S: T^E — * T^E defined as follows: 

• The Liouville section, which is the vertical section corresponding to the 
Liouvillc dilation vector field: 

A(a)=r(a,a) = (a,0,ar)- 

We also mention that the complete lift of a section X G Scc(i?) is the 
section of E characterized by the following properties: 

(i) projects to X, i.e., Tt o X^ = X ot, 

(ii) i-xc-A ^xfJ', 

where by a G C°°{E) we denote the linear function associated to a G Sec(iJ*). 

Given a Lagrangian function L G C°°{E) we define the Cartan 1-section 
0L e Sec((T^£')*) and the Cartan 2-section ujl e Scc{a'^ (T'^E)*) and the 
Lagrangian energy El G C°°{E) as 

eL^S*{dL), ujL^-d0L and El=I:^aL-L. (3.1) 

If (x',?/") are local fibred coordinates on E, (Pq,C'2^) are the corresponding local 
structure functions on E and {Xq,,Vq,} is the corresponding local basis of sections 
of T^E then 

S%a = V„, SVc, = 0, for aU a, (3.2) 
A = 2/"V„, (3.3) 



Br 

El = g^y- - L. (3.5) 

From (jnSI), and it follows that 

isxi^L = -S*{ixi^L), lAi^L ^ -S*{dEL), (3.6) 

for X G Scc{T^E). 

Now, a curve t — > c(t) on i? is a solution of the Euler- Lagrange equations 



for L if 



c is admissible (that is, p{c{t)) = fn{t), where m = t o c) and 
Hc(t),c(t))WL(c(t)) - dEL{c{t)) = 0, for aU t. 
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If c{t) = (a;'(<), y°'{t)) then c is a solution of the Euler-Lagrange equations for L if 
and only if 

Note that if E is the standard Lie algebroid TM then the above equations are 
the classical Euler-Lagrange equations for L : TAI — * R. 

On the other hand, the Lagrangian function L is said to be regular if is a 
symplectic section, that is, if ujl is regular at every point as a bilinear form. In 
such a case, there exists a unique solution verifying 

In addition, using (|3.t)|) . it follows that isTl^l ~ iA^L which implies that is 
a SODE section, that is, 

or alternatively TT{TL{a)) ~ a for all a ^ E. 

Thus, the integral curves of (that is. the integral curves of the vector field 
P^(^l)) arc solutions of the Euler-Lagrange equations for L. F^ is called the 
Euler-Lagrange section associated with L. 



dy^dy 
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From H,S.4|I , we deduce that L is regular if and only if the matrix Wap 
is regular. Moreover, the local expression of F^ is 

where the functions /" satisfy the linear equations 
Examples. 

1. - Real Lie algebras of finite dimension. Let g be a real Lie algebra of 
finite dimension and L : g — > M be a Lagrangian function. Then, the Euler- 
Lagrange equations for L are just the well-known Euler-Poincare equations 
(see, for instance, |36p. 

2. - The tangent bundle. Let L : TM — > M be a standard Lagrangian function 
on the tangent bundle TM of AL Then, the resultant equations are the classical 
Euler-Lagrange equations for L. 

3. - Foliations. If the Lie algebroid is the tangent bundle of a foliation J' on P 
then one recovers the classical formalism of holonomic mechanics. 

4. - Atiyah algebroids. Let tq\G : TQ/G — > M be the Atiyah algebroid 
associated with a principal G-bundle p : Q — > M and L : TQ\G — > M be 
a Lagrangian function. Then, the Euler-Lagrange equations for L are just the 
Lagrange-Poincare equations (see |2()|). 

5. - Action Lie algebroids. Suppose that g is a real Lie algebra of finite di- 
mension and that $ : g x 1/* — > V* is a linear representation of g on V*. If 
L : QxV* — > R is a Lagrangian function on the action Lie algebroid gxV* — > V* 
then the Euler-Lagrange equations for L are just the so-called Euler-Poincare 
equations with advected parameters or the Euler-Poisson-Poincare equa- 
tions (see \20\). 
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3.2. Hamiltonian Mechanics. In this section, we discuss how the Hamiltonian 
formahsm can be developed for systems evolving on Lie algebroids (for more details, 
see 123 EHl)- 

Let T* : E* — > M be the vector bundle projection of the dual bundle E* to E. 
Consider the prolongation T^E* of E over r*, 

T'^E* = {{h,v) eE-KTE* \ p{b) = {TT*){v)} 

= {{a*,b,v) e E* X E xTE* \ t* {a*) ^ T{b), p{b) = {Tt*){v)} . 

The canonical geometrical structures defined on T^E* are the following: 



The Liouville section Qe £ Scc((T^i?*)*) defined by 

QE{a*)ib,v) = a*{b). (3.9) 

3 canonical sym 

by 



• The canonical symplectic section Qe G Sec(A^(T^£'*)*) is defined 



nE = -deE- (3.10) 

where d is the differential on the Lie algcbroid T^E* . 

Take coordinates (a;*,pa) on E* and denote by {^a, T^} the local basis of sections 
r^S*, with 

y„(a*)= (^a*,e„(r*(a*)),pL^) and ^/^(a*) = (^a*, 0, ^) . 
In coordinates the Liouville and canonical symplectic sections are written as 



Qe = PcT and = A + ^P^Cl^^T A 



where {T^^t}} is the dual basis of {^a, J"^}- 

Every function H € C°° {E* ) define a unique section Th of T^E* by the equation 

ith^e = dH, 

and. therefore, a vector field p"^ (Th) = Xh on E* which gives the dynamics. In 
coordinates, 

dH ( . dH 9H\ 



and therefore. 



_ ,dH d f ,dH dH\ d 

Thus, the Hamilton equations are 

dx' _ dH dpa _ , dH dH 

The Poisson bracket { , defined by the canonical symplectic section fi^; on 
E* is the canonical Poisson bracket, which is known to exists on the dual of a Lie 
algebroid [3]. 
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Examples. 

1. - Real Lie algebras of finite dimension. If the Lie algebroid i5 is a real 
Lie algebra of finite dimension then the Hamilton equations are just the well-known 
Lie-Poisson equations (see, for instance, |36p. 

2. - The tangent bundle. If E is the standard Lie algebroid TM and H : 
T*M — > R is a Hamiltonian function then the resultant equations are the classical 
Hamilton equations for H . 

3. - Foliations. If the Lie algebroid is the tangent bundle of a foliation J then 
one recovers the classical formalism of holonomic Hamiltonian mechanics. 

4. - Atiyah algebroids. Let tq\G : TQ/G — > M ^ Q\G he the Atiyah alge- 
broid associated with a principal G-bundlc p : Q — > M and H ; T*Q\G — > M be 
a Hamilton function. Then, the Hamilton equations for H are just the Hamilton- 
Poincare equations (see |26j'). 

5. - Action Lie algebroids. Suppose that g is a real Lie algebra of finite dimen- 
sion, that y is a real vector space of finite dimension and that ^ : qxV* — > V* is 
a linear representation of g on y*. li H : g* x V* — > R is a Hamiltonian function 
on the action Lie algebroid q x V* — > V* then the Hamilton equations for H are 
just the Lie-Poisson equations on the dual of the semidirect product of 
Lie algebras s = q(§)V (see |20|'). 

3.3. The Legendre transformation and the equivalence between the La- 
grangian and Hamiltonian formalisms. Let L : E — > M be a Lagrangian func- 
tion and 9l G Scc((T^£^)*) be the Poincare-Cartan 1-section associated with L. 

Wc introduce the Legendre transformation associated with L as the smooth 
map LegL ■ E — > E* defined by 

LegL{a){b)^j^L{a + tb)l^^, (3.12) 

for a,b G E^, where Ex is the fiber of E over the point x E M. In other words 
LegL{a){b) = 9L{a){z), where 2 is a point in the fiber of E over the point a such 
that Tt(z) = h. 

The map LegL is well-defined and its local expression in fibred coordinates on 
E and E* is 

BT 

LegL{x\y'^) = {x\—). (3.13) 

From this local expression it is easy to prove that the Lagrangian L is regular if 
and only if Leg^ is a local diffeomorphism. 

The Legendre transformation induces a map TLegL '■ T^E — > T^E* defined 

by 

{TLegL){b,Xa) = {b, iTaLegL){Xa)), (3.14) 

for a,b€ E and {a,b,Xa) <E TfE C Er(^a) x Er(a) x TaE, where TLegL ■ TE — > 
TE* is the tangent map of LegL- Note that r* o LegL = t and thus TLegL is 
well-defined. 

If wc consider local coordinates on T^ E (rcsp. T^E*) induced by the local basis 
{Xq, Vq} (resp., {^Q, CP"}) the local expression of TLegL is 

TLegLix\y-;z-,v-) = {x\^;z-,p},z^^^+v^'^^). (3.15) 

The relationship between Lagrangian and Hamiltonian Mechanics is given by 
the following result. 
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Theorem 3.1. |26j The pair [TLegL^LegL) is a morphism between the Lie al- 
gebroids (T^£;, I-, -f , p^) and {T^ E* -f ', P"'*)- Moreover, if Ol and ujl (re- 
spectively, Qe cind ^Ie) CLfe the Poincare-Cartan l-section and 2-section associated 
with L (respectively, the Liouville l-section and the canonical symplectic section on 
T^E*) then 

{TLegL,LegLT{QE) = Ol, (TLegL, LegL)* {ns) = lol- (3.16) 

In addition, in (2^, it is proved that if the Lagrangian L is hyperregular, that 
is, LegL is a global diffeomorphism, then {TLegL, LegL) is a symplectomorphism 
and the Euler-Lagrangc section associated with L and the Hamiltonian section 
Th are (TLegL, Leg L)-"cclatcd, that is, 

Th o LegL = TLegL ° ^l- (3-17) 

Therefore, an admissible curve a{t) on E is a solution of the Euler-Lagrange 
equations if and only if the curve fi{t) = LegL{a,{t)) is a solution of the Hamilton 
equations. 

4. NONHOLONOMIC LAGRANGIAN SYSTEMS ON LlE ALGEBROIDS 

4.1. Constrained Lagrangian systems. In this section, we will discuss La- 
grangian systems on a Lie algebroid r : E — > M subject to nonholonomic con- 
straints. The constraints are real functions on the positions and generalized veloci- 
ties which constrain the motion to some submanifold M of i?. M is the constraint 
submanifold. 

We will assume that the constraints are purely nonholonomic, that is, not all the 
generalized velocities are allowable, although all the positions are permitted. So, 
we will suppose that tt ~ Tjjvt : M — > M is a fibration. 

The constraints arc linear if they arc linear functions on E or, in more geometrical 
terms, if M is a vector subbundle of E over M (Lagrangian systems subject to linear 
constraints were discussed in [r^HS] ). 

In the general case, since tt is a fibration, the prolongation T^M is defined. 
We will denote by r the dimension of the fibers of tt: M — > M, that is r = 
dimM - dimM. 

Now, we define the bundle V — > M of virtual displacements as the subbundle 
of T*E of rank r whose fiber at a point a G M is 

Va = { 6 G Er(a) \ G T„M } . 

In other words, the elements of V are pairs of elements (a, b) ^ E ®m E such that 



^(^(a + tb) 



= 0, 

t=o 



for every local constraint function (/). 

We also define the bundle of constraint forces \I/ by = S'*((T^M)°). Since 
TT is a fibration, the transformation S* : (T^M)*^ — > ^ defines an isomorphism 
between the vector bundles (T^M)" — > M and 'S' — > M. Therefore, the rank of 
^ is s = n — r, where n is the rank of E. 

Next, suppose that L G C°°{E) is a regular Lagrangian function. Then, the pair 
(L, M) is a constrained Lagrangian system. Moreover, assuming the validity 
of a Chetaev's principle in the spirit of that of standard Nonholonomic Mechanics 
(see |2S]), the solutions of the system (L, M) are curves t — > c{t) on E such that: 

- c is admissible (that is, p{c(t,)) = rn{t), where to = r o c), 

- c is contained in M and, 
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- i^^^t)At))'^L{c{t)) - dEL{c{t)) e *(c(i)), for all t. 

If (a;*,y") arc local fibred coordinates on ^^"^ ^hc corresponding 

local structure functions of E and 

are the local equations defining M as a submanifold of E, then { — — X" > 

L aj/" J A=i,...,s 

is a local basis of ^ . Moreover, a curve t — * c{t) = (a;*(t), y"{t)) on _E is a solution 
of the problem if and only if 

0^(x\y") = O, 
where Xa are the Lagrange multipliers to be determined. 

These equations are called the Lagrange- d'Alembert equations for the con- 
strained system (i,M). Note that if E is the Lie algebroid TM, then the above 
equations are just the standard Lagrange-d'Alembert equations for the constrained 
system (L, M). 

Now, we will assume that the solution curves of the problem are the integral 
curves of a section F of T^E — > E. Then, we may reformulate geometrically the 
problem as follows: we look for a section F of T^E — > E such that 

{iruJL - dEL)\M eSeci"^), 

(4.2) 

F|m G Scc(T^M). 
If F is a solution of the above equations then, from (|3.6|l . we have that 

(isr^L - jAWL)|M = 0, 

which implies that F is a SODE section along M, that is, {ST — A)|M = 0. 

4.2. Regularity, projection of the free dynamics and nonholonomic bracket. 

We will discuss next the regularity of the constrained system {L, M) (the con- 
strained system (L,M) is regular ii equations (|4.2|l admit a unique solution F). 

For this purpose, we will introduce two new vector bundles F and T^M over M. 
The fibers of F and T^M at the point a e M are 

Fa = coZ'i^a), 

rjM ^ {ze TJ^M I Tiriz) eVa} = {ze T/M | S{z) e T/M } . 

Then, one may prove the following result. 

Theorem 4.1. |1()| The following properties are equivalent: 

(i) The constrained Lagrangian system (L,M) is regular. 

(ii) r^MnF = {0}. 

(iii) r^Mn(r^M)^ = {o}. 

Here, the orthogonal complement is taken with respect to the symplectic sec- 
tion lul. 

Condition (ii) (or, equivalently, (iii)) in Theorem 14. II is locally equivalent to the 
regularity of the matrix 

V 9?/" dyP ' A,B=i,...,s 
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where {W^^) is the inverse matrix of ( Wap = 



Thus, if L is a Lagrangian function of mechanical type (that is, L[a) = ^S{a,a)^ 
F(r(a)), for aU a^E, with S : ExmE — > R a bundle metric on E and V : M — > 
R a real function on Af ) then the constrained system (L, M) is always regular. 

Now, assume that the constrained Lagrangian system {L, M) is regular. Then 
(ii) in Theorem is equivalent to {T^E)\jy( = T^M ® F and we wiU denote by 
P and Q the complementary projectors defined by this decomposition 

Pa-.T^E ^T^M, Qa-.TfE—^Fa, for all a eM. 

Moreover, we have 

Theorem 4.2. |1(J) Let (L, M) be a regular constrained Lagrangian system and let 
Tl be the solution of the free dynamics, i.e., ivL^L = dE^. Then, the solution of 
the constrained dynamics is the SODE F along M obtained as follows 

r = p(rL|M). 

On the other hand, (3) in Theorem lO is equivalent to {T^E)\j^ = T^M 8 
(T^M)-*" and we will denote by P and Q the corresponding projectors induced by 
this decomposition, that is. 

Pa : TfE — > rjM, Qa = TfE — > (TJM)^, for all a G M. 

Theorem 4.3. ,10' Let (L,M) be a regular constrained Lagrangian system, 
(respectively, T) be the solution of the free (respectively, constrained) dynamics and 
A be the Liouville section ofT^E — > E. Then, T = P{Ti^\'M) if and only if the 
restriction to M of the vector field p^(A) on E is tangent to M. 

Note that if M is a vector subbundlc of E then the vector field p^(A) is tangent 
to M. Therefore, using Theorem 14.31 it follows that 



Corollary 4.4. Under the same hypotheses as in Theorem \4-S\ if M is a vector 
subbundlc of E (that is, the constraints are linear) then F = P(ri|M). 

Next, we will study the conservation of the Lagrangian energy for the constrained 
Lagrangian system (i,!M). 

Since S* : (T^M)" — > is a vector bundle isomorphism, it follows that there 
exists a unique section Q;(ijvt) of (T^M)" — * M such that 

Moreover, we have 

Theorem 4.5 (Conservation of the energy). ^Uj Lf {L,M.) is a regular constrained 
Lagrangian system and F is the solution of the dynamics then iLr(-E'L|M) ~ if 
and only i/a(L jvt)(^l^) ~ 0. Ln particular, if the vector field p^{A) is tangent to 
M then LriEL\M) = 0. 

Now, suppose that / and g are two smooth real functions on JVt and take ar- 
bitrary extensions to E denoted by the same letters. Then, we may define the 
nonholonomic bracket of / and g as follows 

{f,g}r.h=UJL{P{Xf),P{Xg))\^, 

where Xf and Xg are the Hamiltonian sections on T^E associated with / and g, 
respectively. 

The nonholonomic bracket is well-defined and, furthermore, it is not difficult to 
prove the following result. 
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Theorem 4.6 (The nonholonomic bracket). |10| The nonholonomic bracket is an 
almost-Poisson bracket, i.e., it is skew- symmetric and satisfies the Leibniz rule (it 
is a derivation in each argument with respect to the usual product of functions). 
Moreover, if f € C°°(M) is an observable, then the evolution f of f is given by 

f^p^RM) + {f,EL\MUh, 

where Rl is the section o/T^M — > M defined by Rl = P{Vl\'M) - P{Tl\'M). In 
particular, if the vector field {l!^) is tangent to M then 

/ = {/,i?i|M}„,,. 

4.3. Reduction. Next, we will discuss a reduction process and its relation with Lie 
algebroid epimorphisms. These results will be also valid for Lie algebroids without 
nonholonomic constraints, simply taking M = i? and M' = E' in the sequel. 

Let (L, M) be a regular constrained Lagrangian system on a Lie algebroid r : 
E — > M and (L', M') be another constrained Lagrangian system on a second Lie 
algebroid r' : E' — * AI' . Suppose also that we have a fiberwise surjective morphism 
of Lie algebroids (f> : E — > E' over a surjective submersion cj) : M — > M' such 
that: 

i) L = L'o$, 

ii) $|jv( : M — > M' is a surjective submersion and 
zm) $(V,)=V:,(^),forallaGM. 

Note that if M and M' are vector subbundles of E and E' , respectively, then 
conditions z), ii) and Hi) hold if and only if 

L = L' o $ and $(M) = M'. 

In the general case, one may introduce the map T*(f> : T-^M — > M' given by 

(T*$)(6, v) = ($(6), for (6, v) e T^M, 

and we have that T*$ is a Lie algebroid epimorphism over <I>. In addition, the 
following results hold 

Theorem 4.7 (Reduction of the constrained dynamics). ^0] Let (L,M) be a reg- 
ular constrained Lagrangian system on a Lie algebroid E and {L','M') be a con- 
strained Lagrangian system on a second Lie algebroid E' . Assume that we have a 

fiberwise surjective morphism of Lie algebroids $ : E > E' over <p '■ M — > M' 

such that conditions i), ii) and Hi) hold. Then: 

(i) The constrained Lagrangian system (L',M') is regular. 

(ii) If T ( respectively, T' ) is the constrained dynamics for L ( respectively, for 

L') then r*$or = r'o$. 

(iii) If t — > c{t) is a solution of Lagrange-d'Alembert differential equations for 
L then <E>(c(t)) is a solution of Lagrange-d'Alembert differential equations 
for V. 

Theorem 4.8 (Reduction of the nonholonomic bracket). JO] Under the same 
hypotheses as in Theorem \4. T[ we have that 

{/'o$,.g'o$}„,, = {/',g'};,o<i> 

for f',g' G C°°(M'), where {■,-}nh (respectively, {-j-Ynh) nonholonomic 
bracket for the constrained system (L,M) (respectively, (L',M')j. In other words, 
4> : M — > M' is an almost-Poisson morphism when on M and M' we consider the 
almost-Poisson structures defined by the corresponding nonholonomic brackets. 
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Reduction by symmetries. Let (/) : Q — > M be a principal G-bundle and r : 
E — > Q be a Lie algebroid over Q. In addition, assume that we have an action of 
G on E such that the quotient vector bundle E/G is defined and the set Sec{E)'~^ 
of equivariant sections of i? is a Lie subalgebra of Sec{E). Then, E' = E/G 
has a canonical Lie algebroid structure over M such that the canonical projection 
$ : E — > E' is a fiberwise bijective Lie algebroid morphism over (j) (see Uni). 

Next, suppose that (L,M) is a G-invariant regular constrained Lagrangian sys- 
tem, that is, the Lagrangian function L and the constraint submanifold M are 
G-invariant. Assume also that M is closed. Then, one may define a Lagrangian 
function L' : E' > R on E' such that 

Moreover, G acts on M and the set of orbits M' = M/G of this action is a quotient 
manifold, that is, M' is a smooth manifold and the canonical projection $|jvt : 
M — > M' = M/G is a submersion. Thus, one may consider the constrained 
Lagrangian system (i',M') on E' . 

Since the orbits of the action of G on i? are the fibers of $ and M is G-invariant, 
we deduce that 

Va($) C TaM, for aU a £ M, 
being the vertical bundle of This imphes that $|v„ : — > '^*(a) ^ 
linear isomorphism, for all a 6 M. 

Therefore, from Theorem 14. 71 we conclude that the constrained Lagrangian sys- 
tem (i',M') is regular and that 

r*<I> o F = F' o $, 

where F (resp., F') is the constrained dynamics for L (resp., L'). In addition, using 
Theorem l4.8l we obtain that $ : M — > M' is an almost-Poisson morphism when on 
M and M' wc consider the almost-Poisson structures induced by the corresponding 
nonholonomic brackets. 

4.4. Example: a rolling ball on a rotating table. We apply the results in this 
section to the case of a ball rolling without sliding on a rotating table with constant 
angular velocity [Tl HI [TT) ! 1301 148j . A (homogeneous) sphere of radius r > 0, unit 
mass TO = 1 and inertia about any axis, rolls without sliding on a horizontal 
table which rotates with constant angular velocity about a vertical axis through 
one of its points. Apart from the constant gravitational force, no other external 
forces are assumed to act on the sphere. 

Choose a Cartesian reference frame with origin at the center of rotation of the 
table and z-axis along the rotation axis. Let (x, y) denote the position of the point 
of contact of the sphere with the table. The configuration space for the sphere on 
the table is Q = x SO{3), where 50(3) may be parameterized by the Eulerian 
angles 9, ip and ^j. The kinetic energy of the sphere is then given by 

T = ^(.T^ +f + k\d^+ + 2ipijj cos 6*)). 

With the potential energy being constant, we may put V ~ 0. Thus, the Lagrangian 
function L is T and the constraint equations are 

X — r 8 sinip + rep sin 9 cos "ip = — Hy, 
y + r6 cos Ip + np sin 9 simjj ~ fix. 

Since the Lagrangian function is of mechanical type, the constrained system is 
regular. Note that the constraints are not linear and that the restriction to the 
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constraint submanifold M of the Liouville vector field on TQ is not tangent to M. 
Indeed, the constraints are linear if and only if f2 = 0. 

Next, following we will consider local coordinates (x, y, ^, -0; 7ri)i=i_..._5 

on TQ = TM? x T(S'0(3)), where 

x = x, y = y, = 0, 'f = (p, '0 = 1/', 
m = rx + k'^q2, 7r2=ry-fc^gi, tts = fc^ga, 

^2 

^4 = (fc2 ^ ,.2) - '■'?2 + ^y)^ ^5 = (^2 ^ p2) + '''^1 - ^^2^), 

(gi, ^2, 93) are the quasi-coordinates defined by 

qi = UJx, 92 = (^y, 93 = 
and w^,, and are the components of the angular velocity of the sphere. 
Then, the constrained dynamics is the SODE F along M defined by 

r = (pri|M) = (x^ + y± + e-t + ^^ + i^^)\M 

ax ay 80 dip dib n\ 
,. d .8 .8 .8 . 9 

ax ay aqi aq2 8q3 

On the other hand, when constructing the nonholonomic bracket on M, we find 
that the only non-zero fundamental brackets are 

{x,TTi}nh = r, {y,'^2}nh = r, 

{91, 7r2}n/i = — 1, {92, 7I'l}n/i = 1, {93 , ""S } n/i = 1 , 

{T^l,T^2}nh = TTg, {n2,Tr3}nh = ^^2^y2) ^l + (fc2^_^^2) ^^' ^^'"^^ 

in which the "appropriate operational" meaning has to be attached to the quasi- 
coordinates qi. 

Thus, we have that 

/ = RlU) + {f,L}^h, for / e C-(M), 
where is the vector field on M given by 

k^n ril 8 8 

^ ~ ^ {k^ + r^y'' dy ' ^ dx' ^ (k^ + r^y"^ 8qi ^ y dq2 

-fa;(^3 - fc'^^)^ + yi^3 - k'n)— - eiir.x + n2y)^^))\M. 

Note that i?L if and only if f7 = 0. 

Now, it is clear that Q = M.^ x 5*0(3) is the total space of a trivial principal S0{3)- 
bundle over and the bundle projection cj) : Q — > M = is just the canonical 
projection on the first factor. Therefore, we may consider the corresponding Atiyah 
algebroid E' = TQ/SO{3) over M = R^. 

One may prove that E' is isomorphic to the real vector bundle TM^ x M.^ — > 
in such a way that the anchor map p' : E' = TR^ x — > TR^ is just the canonical 
projection on the first factor. Moreover, one may choose a global basis {e^}i=i_..._5 
of Sec(_B') and the only non-zero fundamental Lie brackets are 

[64,41' = 4, 14,64]' = 4, 14,41' = 4- 

We have that the Lagrangian function L = T and the constraint submanifold 
M are S'0(3)-invariant. Consequently, L induces a Lagrangian function L' on 
E' = TQ/SO{3) and the set of orbits M' = M/SO{'S) is a submanifold of E' = 



A SURVEY OF MECHANICS ON LIE ALGEBROIDS AND GROUPOIDS 



23 



TQ / SO{3) in such a way that the canonical projection $|M : M — > M' = 
M/50(3) is a surjective submersion. 

Under the identification between E' = TQ/SO{3) and TM^ x R3, L' is given by 

1 P 

L'{x,y,x,y;u;i,uJ2,uj3) = ^(i^ + y^) + y(t^? + t^2 +^3)' 

where (x,y,x,if) and (cj 1,0^2,^^3) are the standard coordinates on TM? and K"^, 
respectively. Moreover, the equations defining M' as a submanifold of TR^ x 'E? 
are 

i — raj2 + = 0, y + rtJi — f2a: = 0. 

So, we have the constrained Lagrangian system (L',M') on the Atiyah algebroid 
E' = TQ/SO{3) = TR2 X R^. Note that the constraints are not linear and that if 
A' is the Liouville section of the prolongation E' then the restriction to M' of 
the vector field {p'Y (A') is not tangent to M'. 
Now, if we put 

x' = x, y' = y, 

tt[ = rx + k'^uj2, TTj = rij — k'^LOi, TTg ~ fc^W3, 

^2 ^2 

< = (fc2 + ^2) - ''^2 + TT,^ = j^^—-^{y + TLJl - ftx) , 

then (x', tt'j^, TTj, TTg, 7r4, TTg) is a system of global coordinates on TQ/ SO{3) = 
TM^ X M'^. In these coordinates the equations defining the submanifold M' are 
TT^ = and TTg = and the canonical projection $ : TQ — > TQ / SO{3) is given by 

^{x,y,0,<f,{i;;TTi,7T2,Tr3,Tr4,7T5) = (i,y;7ri,7r2,7r3,7r4,7r5). (4.5) 

Thus, if r' is the constrained dynamics for the system (i',M'), it follows that 
(see (1331) 

On the other hand, from (|4.4|l . (|4.5() and Theorem 14.81 we deduce that the only 
non-zero fundamental nonholonomic brackets for the system (L',!M') are 

I'ri,7r2}„,^ ^TTg, {'^1.T^3]nh ^ (p^^2)^l + (A;2+^2)2^ ' 

i'^3-'^ll„/.- (p_^^2)'^2 (fc2+^2)^- 

Therefore, we have that 

/' = {(^Y\rM) + for /' e C-(M'), 

where (p')"^ (^L') is the vector field on M' given by 

4.5. Hamiltonian formalism. Let (£,M) be a constrained Lagrangian system 
on a Lie algebroid E and assume that the Lagrangian function L is hyperregular. 
Then, since the Legi^ is a diffeomorphism then, it is clear that one may develop 
a Hamiltonian formalism which is equivalent, via the Legendre transformation, to 
the Lagrangian formalism. 
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5. Mechanical control systems on Lie algebroids 

5.1. General control systems on Lie algebroids. Consider a Lie algebroid 
r: E — > M, with anchor map p: E — > TM . Let {<t, ryi, . . . , ry^} be sections of E. 
A control problem on the Lie algebroid t: E — > M with drift section a 
and input sections 771, . . . , 77^ is defined by the foUowing equation on M, 

fc 

m{t) = p(c7{m{t)) + u,(t)77,(m(t))) , (5.1) 

i=l 

where u = {ui, . . . , Uk) S J7, and U is an open set of R'^ containing 0. The function 
1 1— > u{t) = (Mi(t), . . . , Uk{t)) belongs to a certain class of functions of time, denoted 
by It, called the set of admissible controls. For our purposes, we may restrict the 
admissible controls to be the piecewise constant functions with values in U . Notice 
that the trajectories of the control system are admissible curves of the Lie algebroid, 
and therefore they must lie on a leaf of E. It follows that if E is not transitive, 
then there are points that cannot be connected by solutions of any control system 
defined on such a Lie algebroid. In particular, the system (|5.1|) cannot be locally 
accessible at points m e M where p is not surjective. Since the emphasis here is 
put on the controllability analysis, without loss of generality we will restrict our 
attention to locally transitive Lie algebroids. 

Denoting by / = p{a) and gi — p{i]i), i £ 1, . . . , fc, we can rewrite the system ()5.1|l 

as 

k 

m{t) = f{m{t)) + Y,Mt)9^{m{t)), (5-2) 
1=1 

which is a standard nonlinear control system on AI afhne in the inputs |47| . Here 
we make use of the additional geometric structure provided by the Lie algebroid 
in order to carry over the analysis of the controllability properties of the control 
system ()5.1|l . We refer to 021 for a comprehensive discussion of the notions of 
reachable sets, accessibility algebra and computable accessibility tests. 

Definition 5.1. The accessibility algebra T) of the control system (|5.1(l in 
the Lie algebroid is the smallest subalgebra of Scc{E) that contains the sections 
{cr,77l, . . .,f]k}- 

Using the Jacobi identity, one can deduce that any element of accessibility alge- 
bra 2) is a linear combination of repeated Lie brackets of sections of the form 

I0,[0-i,I...,IC2,Cil---lll, 

where Q S {cr, rji, . . . ,rik}, I < i < I and I £ N. 

Definition 5.2. The accessibility subbundle in the Lie algebroid, denoted 
by Lie({(T, ?7i, . . . , 77fc}), is the vector subbundle of E generated by the accessibility 
algebra D , 

Lie({cr, 7]i, . . . , r^fc}) = span {C{m) \ C section of E in ©} , m G M. 

If the dimension of Lie({CT, 771, ... , 77^}) is constant, then Lie({(T, 771, ... , 77^}) is 
the smallest Lie subalgebroid of E that has {cr, 771, . . . , 77^} as sections. 

5.2. Mechanical control systems. Let r: E — > M be a Lie algebroid, let V be 
a connection on E, and let {7/, 7/1, . . . , 7;/c} be sections of E. A mechanical control 
system on the Lie algebroid t : E — * M is defined by the following equation 

fc 

^a(t)a{t) + r]{m{t)) = ^ u,{t)7j,{m{t)). (5.3) 

1=1 
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We will often refer to 77 as the potential energy term in equations H5.3|l . Associated 
with this equation, there is always a control system on the Lie algebroid T^E — > E 
given by 

k 

a{t) = ((Fv - ^")(a(i)) + ^ u.iMHt))) , (5-4) 

i=l 

where 77^ (resp. -q]^) denotes the vertical hft of 77 (resp. -qi) and Fy is the SODE 
section associated with V. Fy is locally given by (see 

Fv=y"X„~i(F^^ + F^^) yVV„. 

There are two distinguished families within the class of mechanical control sys- 
tems. We introduce them next. 

Mechanical control systems. Consider a Lagrangian system with L: E — > M 
of the form 

L{a) = -'5{a,a) -V oT{a), aeE, 

with S- E Xm E — > M a bundle metric on E and V a function on M. This 
Lagrangian function gives rise to the Euler-Lagrange equations as explained in 
Section rm 

Consider now the situation when the Lagrangian system is subject to some ex- 
ternal forces, represented by a collection {61, . . . , 9k} of sections of E*. Denote by 
{■qi, . . . ,r]k} the input sections of E determined by the control forces . . . ,6k} via 
the metric, i.e., Oi{X) = S(?7i,X) for all X e Sec(£'). If Fys denotes the SODE sec- 
tion associated with the Levi-Civita connection V^, the controlled Euler-Lagrange 
equations can be written as 

k 

a{t) - p-(Fvs(a(i)) - (gradg Vr{a{t)) + MM Ht))) . (5.5) 

i=l 

Here gradg ^ is the section of E characterized by S(gradgy,X) = dV{X) for 
X G Sec{E). Note that system (|5.5|l is a control problem on the Lie algebroid 
T^E — > E as defined in Section Locally, the equations can be written as 

i=l 

where (Sa/3) are the components of the metric S and (S"'') is the inverse matrix of 

Alternatively, one can describe the dynamical behavior of the mechanical control 
system by means of an equation on E via the covariant derivative. An admissible 
curve a : i I— > a(t) is a solution of the system (|5.5() if and only if 

k 

VS(^)a(t) + gradg V{m{t)) = ^ ?/,(t)r,,(m(i)). (5.6) 

i=l 

This equation corresponds to a mechanical control system (|5.3|l with connection 
V = and sections {gradg V,r]i, . . . , 77^}. 

Mechanical control systems v^fith constraints. Assume a mechanical control 
system with data (S, V, {61, . . . , 9k}) is subject to the constraints determined by a 
subbundlc D of E. Consider the orthogonal decomposition E = D Q) an the 
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associated orthogonal projectors P : E — *■ D, Q: E — > D . Then, one can write 
the controlled Lagrange-d'Alembert equations as 

k 

^(Vf(t)a(t)) + ^(gradg V{m{t))) = Y,u,{t)P{r^,{m{t))), Q{a) = 0. 

1=1 

In terms of the constrained connection Vo-?/ = ^(V^?/) + V^(Q7y), with <T,ri Q 
Sec{E), the controlled equations can be rewritten as Va(j)a(i)+P(gradg V{m{t))) = 

Si=i ''^ii't) P iVii^T-i't))) J Q{o) ~ 0. Since the forcing terms coming from the poten- 
tial and the inputs belong to D, the solutions of the total controlled dynamics 
initially belonging to D also remain in D. As a consequence, an admissible curve 
a :t ^ a{t) is a solution of the system (|5.8() if and only if 

k 

Va(t)a{t) + P(gradg V{m{t))) = ^ u,{t)P{r^,{m{t))), ao £ D. (5.7) 

i=l 

This equation corresponds to a mechanical control system (|5.3|l with connection 
V = V and sections {P(gradg V), P{f]i), . . . , -P(?7fc)}. 

Note that one can write the controlled dynamics as a control system on the Lie 
algebroid T^E — > E, 

k 

a{t) = (r^{a{t)) - P(gradg Vr{a{t)) + ^ u,{t)P{rj,r {a{t))) . (5.8) 

i=l 

The coordinate expression of these equations is greatly simplified if we take a basis 
{cq} = {ea,eA} of E adapted to the orthogonal decomposition E ^ D ® D^, 
i.e., D = spanjca}, 'D^ = spanje^}. Denoting by (?/") = {y'^,y^) the induced 
coordinates, the constraint equations Q{a) = just read y"^ = 0. The controlled 
equations (|5.7|l are then 

= Ply", 

i=l 

where S^^ = Fj!^ + Fg^ are the components of the symmetric product. 

5.3. Accessibility and controllability notions. Here we introduce the notions 
of accessibility and controllability that are specialized to mechanical control systems 
on Lie algebroids. Let m G M and consider a neighborhood V of m in M. Define 
the set of reachable points in the base manifold M starting from m as 

RM(m,r) = {m' e M \3ueU defined on [0,T] such that the evolution of 

for a(0) = Om satisfies r(a(t)) e V, t £ [0,T] and T(a(T)) = m'} . 

Alternatively, one may write IR](f(m,r) = '''(R^ ^^■'(0„i, T)). Denote 

Rl,im,<T)= [j im,t). 

t<T 

Definition 5.3. The system H5.4|l is locally base accessible from m (respec- 
tively, locally base controllable from m) i/ R]^^(?t!,, < T) contains a non-empty 
open set of AI (respectively, M.^{m,< T) contains a non-empty open set of M to 
which m belongs) for all neighborhoods V of m and allT > 0. If this holds for any 
m G M, then the system is called locally base accessible (respectively, locally 
base controllable). 
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In addition to the notions of base accessibility and base controllability, we shall 
also consider full-state accessibility and controllability starting from points of the 
form 0„i G E, m E M (note that full-state is meant here with regards to E, not to 
TM). 

Definition 5.4. The system (15.4(1 is locally accessible from m at zero (re- 
spectively, locally controllable from m at zero) «/M^(Om,< T) contains a 
non-empty open set of E (respectively, M^(0„i,< T) contains a non-empty open 
set of E to which Om belongs) for all neighborhoods W of 0„i in E and all T > 0. 
If this holds for any m £ M, then the system is called locally accessible at zero 
(respectively, locally controllable at zero). 

The relevance of the above definitions stems from the fact that, frequently, one 
needs to control a system by starting at rest. Nevertheless it is important to notice 
that not every equilibrium point at m corresponds to the point Om- Finally, we also 
introduce the notion of accessibility and controllability with regards to a manifold. 

Definition 5.5. Let tp: M — > N be an open mapping. The system 1(5. 4() is lo- 
cally base accessible from m with regards to N (respectively, locally base 

controllable from m with regards to N) if ip{M.\j{m, < T)) contains a non-empty 
open set of N (respectively, '0(R]^^(m, < T)) contains a non-empty open set of N 
to which ip{m) belongs) for all neighborhoods V of m and all T > 0. // this holds 
for any m S M , then the system is called locally base accessible with regards 
to N (respectively, locally base controllable with regards to N). 

Note that base accessibility and controllability with regards to M with idj\/ : 
M — > AI corresponds to the notions of base accessibility and controllability (cf. 
Definition 1^21). Moreover, if the system is base accessible, then it is base accessible 
with regards to TV. The analogous implication for base controllability also holds 
true. 

5.4. Tiie structure of tiie control Lie algebra. The aim of this section is to 
show that the analysis of the structure of the control Lie algebra of affine connection 
control systems carried out in can be further extended to control systems 
defined on a Lie algebroid. The enabling technical notion exploited here is that of 
homogeneity. 

Let i? be a Lie bracket formed with sections of the family X = {Fy, rf( , . . . ,ri^ ^rf^}. 
The degree of B is the number of occurrences of all its factors, and is therefore 
given by 5{B) 8a{B) + S,{B) + ■■■ + 8k{B), where 5o{B), S,{B), ie{l,...,k}, 
and 6k+iiB) correspond, respectively, to the number of times that Fy, rj^" , i £ 
{1, . . . , k}, and r/^ appear in B. For each I, consider the following sets 

Br'(X) = {B bracket inX | 6{B) = /}, Br;(X) = {B bracket inX \ B e Ti}, 

where denotes the set of homogeneous sections of TE of degree I. The notion 
of primitive bracket will also be useful. Given a bracket B in X, it is clear that 
we can write B = [Bi,B2], with Bi brackets in X. In turn, we can also write 
Ba — [Bai,Ba2] for a ~ 1,2, and continue these decompositions until we end up 
with elements belonging to X. The collection of brackets i?2, -Bii, -B12, • ■ • are 
called the components of B. The components of B which do not admit further 
decompositions are called irreducible. A bracket B is called primitive if all of its 
components are brackets in Br_i(X) U Bro(X) U {Fy}. 

Consider the set X' = {Fy — ,ri( , . . . ,77^}. Clearly, the elements in Lie(X') 
are linear combinations of the elements in Lie(X). In fact, for each bracket B' 
of elements in X', let us define the subset S{B') C Br(X) formed by all possible 
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brackets B G Br(X) obtained by replacing each occurrence of Fy — in B' by 
either Fy or if \ Then, one can prove by induction (cf. [5H|) that 

J2 i~l)^''^''-^^B. (5.9) 

BeS(B') 

Reciprocally, given an clement B £ Br(X), one can determine the bracket B' of 
elements in X' such that B G S{B') simply by substituting each occurrence of Fy 
or T]^ in B by Fy — ry^. We denote this operation by pseudoinv(i3) = B'. For each 
fc G N, define the following families of sections in E, 

{a G Sec{E) \ a"' = B", B" = ^ {-lY'^+^iB) b , B G Br^^^-^CX) primitive}, 

BeS(pseudoinv(S)) 
nBr_i(X)nBro(X) 

{a G Sec(£:) | a = = ^ 5 ^ primitive}. 

BeS(pseudoinv(B)) 
nBr_i(X)nBro(X) 

Consider 

ehor(??;'7i, ■ • • ,?/fc) = UfcgNe|j^l(?7;'7i: ■ • ■ ,'7fc), 

and denote by Cver(?7; • ■ • , Vk) and Chor('7; ■ • ■ , ?7fc), respectively, the subbun- 
dles of the Lie algebroid E generated by the latter families. 

Taking into account the previous discussion, we are now ready to compute Lie({Fy- 
r]^ ,r]i , . . . ,r]^}) fov a. mechanical control system defined on a Lie algebroid. 

Proposition 5.6. (^) Let m G M. Then, 

Lii({Fy - 7?^77r, . . . ,77^) n Vero„(ri?) = C^^Mvu ■ ■ • ,^fc)M", 
Lie({Fy - if , 111', ■ ■ ■ ,?7fe}) nHor„(T£') = Chor(?7;??i, • ■ ■,Vk){ni). 

Remark 5.7. In the absence of potential terms, i.e., 77 = 0, one has that 



Cvcr(0;?7i, . . . ,77fc) = Sym({?7i, . . .,77^}), Chor(0;ryi, . . . ,77^) = Lie(Sym({77i, . . .,77^-})), 



where Sym({7;i, . . . ,77^}) denotes the distribution obtained by closing (the distri- 
bution defined by) {r/i, ...,77^} under the symmetric product associated with V. It 
is worth noticing that, in this case, Cvcr(0; 771, ... , 77^) C Choi (0; rji, . . . , 77^). This is 
not true in general. o 

5.5. Accessibility and controllability tests. In this section we merge the no- 
tions introduced in Section 15^ with the results obtained in Section to give tests 
for accessibility and controllability. 

Proposition 5.8. Let m G M and assume the Lie algebroid E is locally 
transitive at m. Then the mechanical control system ()5.4|) is 



• locally base accessible from m if ChoriVi Vii ■ ■ ■ 1 '?fe)('^) + kerp = Em, 

• locally accessible from m at zero if Chor('7;'7ij • ■ • i^fe)("^) + ker p = E„ 



and Cvor(f?;?7i, ■ • ■,'nk)im) = E„ 



In order to state controllability tests, we need to introduce the notions of good 
and bad symmetric products. We say that a symmetric product P in the sections 
{77, 771, . . . , 77fc} is bad if the number of occurrences of each 77^ in P is even. Otherwise, 
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P is good. Accordingly, {r/i : rji) is bad and {{rj : ly) : (rji : 77,)) is good. The 
following theorem gives sufficient conditions for local controllability. 

Proposition 5.9. |13| Let m G M . The mechanical control system H5.4|) is 

• locally base controllable from m if it is locally base accessible from m and 
every bad symmetric product in {77,771, . . . ,77fe} evaluated at m can be put 
as an M.-linear combination of good symmetric products of lower degree and 
elements of ker p, 

• locally controllable from m at zero if it is locally accessible from m at zero 
and every bad symmetric product in {77, 771, . . . , 77^} evaluated at m can be 
put as an W-linear combination of good symmetric products of lower degree. 

The corresponding tests for base accessibility and controllability with regards to 
a manifold can be proved in a similar way. 

Proposition 5.10. [T3| Letip: M — > N be an open map. Let m g AI and assume 
')/'*(p(£'„i)) = T^(jn)N . Then the mechanical control system (|5.4() is 

• locally base accessible from m with regards to N if Chor('7; '7i, ■ • ■ , '7fc)('7i) + 
p~i(kerV'*) = E„i, 

• locally base controllable from m with regards to N if the system is locally 
base accessible from m with regards to N and every bad symmetric product 
in {rj, r]i, . . . ,rik} evaluated at m can be put as an W-linear combination of 
good symmetric products of lower degree and elements o/p^^(ker7/'*). 

6. Discrete Mechanics on Lie groupoids 

In this section, we discuss discrete Lagrangian Mechanics on a Lie groupoid 
G ^ M. Instead of the usual Euler-Lagrange equations (|3.7() for a Lie algebroid 
T : E — > M equipped with a Lagrange function L : E — > R, we obtain a 
set of difference equations called Discrete Euler-Lagrange equations for a discrete 
Lagrangian Ld : G — > R [33] . When the Lie algebroid is precisely E = Eq and Ld 
is a suitable approximation of the continuous Lagrangian L : Eq — > R, then we 
will obtain a geometric integrator for the Euler-Lagrange equations. In the next 
subsections we will carefully analyze this construction and its geometric properties. 

G 

6.1. Lie algebroid structure on the vector bundle tt"^ : Eq x Eq — > G. Let 

G =1 M be a Lie groupoid with structural maps 

a,(3:G — > M, e : M — > G, i : G — > G, 777 : G2 — > G. 

G 

We know that Eq x Eg ^ Eg is a Lie groupoid. The following theorem shows that 

G 

the vector bundle tt"^ : Eq x Eg = V (3 ©g Va — > G is equipped with a natural 
structure of Lie algebroid. 

G 

Theorem 6.1 (See Theorem 3.3 in |2S|)- The vector bundle tt"^ : Eg x Eg = 
V (3 ®G Vol — > G admits a Lie algebroid structure, where the anchor map is given 
by 

pEG'i.Ea^^Xg,Yg)^Xg+Yg, for{Xg,Yg)^VgP®Vga, (6.l) 

G 

and the Lie bracket |-, ■JEg'xEg ^/jg gpace Sec(7r'^) is characterized by the following 
relation 

l(X, Y), (X', Y')fo'iEa ^ (-lx,xi iY,Y'i} , (6.2) 
for X,Y,X',Y' e Sec(r). 
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We also remark that, if we denote by (G, a) the fibration a: G — > M and 
by {G,P) the fibration /3: G — > M, then it is not difficult to prove that the Lie 
algebroid prolongations T^^ (G, a) — > G and T^"^ (G, /3) — > G are isomorphic, as 
Lie algebroids, and both are isomorphic to the Lie algebroid V(3 (Bg Vo: — > G and 

G 

hence to Eg x Eg — > G (see [23] for the details). 

G 

The following diagram shows both structures of Eg x Eg, 

^Eg 



Eg X Eg 



a 




TG 




TP 








r 





TM 



=^ M 

where the vertical maps are morphisms of Lie groupoids and the horizontal maps 
are morphisms of Lie algebroids. 

Given a section X of Eg — > M, we define the sections X'-^'°\ (the /3 and 

G 

a- fifts) and X^^'^^ (the complete fift) of X to t:^ : Eg x Eg — > G as follows: 
X(i'0)(5) = (X(5),Og), X("'i)(g) = (Og,^(g)) and X^^^'^g) ^ {~X {g),X {g)) 
We can easily see that 



,0) 



and =0 (6.3) 



and, as a consequence, 



(6.4) 



6.2. Discrete Variational Mechanics on Lie groupoids. Discrete Lagrangian 
systems on Lie groupoids have a variational origin, as we explain next. A discrete 
Lagrangian system consists of a Lie groupoid G ^ AI (the discrete space) and 
a discrete Lagrangian Ld : G — > M. 

Discrete Euler-Lagrange equations. For g E G fixed, we consider the set of 
admissible sequences: 

{ (51, . . . , g^v) e G^ I {gk, 9k+i) e G2 for fc = 1, . . . , iV - 1 and gi . . . g„ = g } . 



It is easy to show that we may identify the tangent space to with 



4si 



g^)^^ ={{vu...,VN-i) I Vk e {Eg)x, andxfc =/3(gfc),l < fc < - 1 } , 



An element of Tj^^ gN)^g called an infinitesimal variation. Now, we define 



the discrete action sum associated to the discrete Lagrangian Ld ; G 



by 



N 



S>Ldiigi 



,9N 



k=l 
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Hamilton's principle requires that this discrete action sum be stationary with re- 
spect to all the infinitesimal variations. This requirement gives the following alter- 
native expressions for the discrete Euler-Lagrange equations (see |H5jl: 

X{gk){Ld)-X{gk+i)iLd)^0, (6.5) 

or 

(did,X(°'^))(5fc) - (dLrf,X(i'«))(gfc+i) = 0, 

for all sections X oi t : Eq — > M. Here, d denotes the differential of the Lie 

G 

algebroid tt"^ : Eq x Eq = V(3 ®c Vet — > G. Alternatively, we may rewrite the 
Discrete Euler-Lagrange equations as 

d° [Ldolg^ +Ldorg^^^ oi] (e(a;fc))|(^^j^^ =0, 

where P{gk) = ce{9k+i) = Xk, and where d° denotes the standard differential on G, 
that is, the differential of the Lie algebroid : TG — > G. 

Thus, we may define the discrete Euler-Lagrange operator: 

DbelLci ■ G2 — > Eq , 
where Eq is the dual of Eq- This operator is given by 

DDEhLdig, h) = d° [Ld ol +LdO rh o i] (e(x))| 

with /3(.g) = a{h) = x. 

G 

Discrete Poincare-Cartan sections. Consider the Lie algebroid tt"^ : Eq x Eq = 
VP (Bg Vet — > G, and define the Poincare-Cartan 1-sections 0£^,0j^ G 
Sec((7r'^)*) as follows 

Ql^{g){Xg,Yg) = -Xg{Ld), et^{g){Xg,Yg) = Yg{Ld), (6.6) 
for each g eG and {Xg,Yg) G VgP Vga. 

Since dLd = Qt — ©7 and so, using = 0, it follows that dOt = ^07 • This 
means that there exists a unique 2-section Hl^ ~ ^'^©Ld ~ ~'^'^Ld' that will be 
called the Poincare- Cartan 2-section. This 2-section will be important to study 
the symplectic character of the discrete Euler-Lagrange equations. 

If {Xi} is a local basis of Sec(T) then {xf"'^\ xf'^'^} is a local basis of Sec(7r'^). 
Moreover, if {{X^)^'^'°\ (X')(0'i'} is the dual basis of '^'l, it follows that 

eZd = -x,{Ld){x^r'''\ e+ = x,{Ld){x'^r''\ 

Hl, = -X,(X,Ld)(X^)(i^o) A (X^)(0'i). 

Discrete Lagrangian evolution operator. Let ^ : G — > G be a smooth map 
such that: 

- graph(^) C G2, that is, {g,£,{g)) € G2, for all g G G (<^ is a second order 
operator). 

- {g,£,{g)) is a solution of the discrete Euler-Lagrange equations, for all g G 
G, that is, (L»QELLd)(5,^(.g)) = 0, for aU g G G. 

In such case 

X{g){Ld)-X{mKLd)^0 (6.7) 

for every section X of Eq and every g ^ G. The map ^ : G — > G is called a 
discrete flow or a discrete Lagrangian evolution operator for Ld- 
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Now, let ^ : G — > G be a second order operator. Then, the prolongation 
: V(3(Bg Vet — > V(3(Bg Vet of ^ is the Lie algebroid morphism over ^ : G — > G 
defined as follows (see 

TgaXg,Yg) = ((T,(r,^(,)Oz))(yg),(r,0(^s) + M(>^s)-r,(r,^(,)Oz)(rg)), (6.8) 

for all {Xg,Yg) € Vgf3®Vga. Moreover, from ifTTji . and we obtain that 

Tg^Xig), Yig)) = (-?($(<?)), (T,C)(X(g) + F (5)) + ¥{^9))), (6.9) 

for all X, Y sections of Eg ■ 

Using 1)6.8(1 . one may prove that (see 

(i) The map ^ is a discrete Lagrangian evolution operator for Ld if and only 

if(n,0*ez^ = e+^. 

(ii) The map ^ is a discrete Lagrangian evolution operator for Ld if and only 
if iT^,0*Ql,-Q-L,=dLd. 

(iii) If ^ is discrete Lagrangian evolution operator then {T^,£^)*flL^ ~ ^II^- 

Discrete Legendre transformations. Given a Lagrangian Ld : G — > R we 
define the discrete Legendre transformations ¥~Ld : G — > Eq and F+L^ : 
G — > by 

(F"Ld)(ft.)(w£(a(,,))) = -We(a(h))(irfOniO«), for We(a(h)) G (£^G)a(/i), 

(F+L)(.g)(i;,(^(g))) = W£(/3(g)) (L<i o Ig), for W£(/3(g)) £ (£^G)/3(ff)- 

Now, we introduce the prolongations TW^Ld : Eq x Eq = Vj3 ©g Va — > T^^ Eq 

and TF+Lrf :EgxEg = VP ©g — > T^g£:^ by 

r^F-Lrf(X;„y;,) = (T^(i or,,-i)(X;,),(T,,F-L)(X,0 + (T,,F-L)(n)), 

7^,F+Lrf(X;„y;,) = inih-iiYh), {Th¥+L){Xh) + (TftF+L)(n)), 

for all /i G G and (X/i, Y/j) G VhP © V;ia. We observe that the discrete Poincare- 
Cartan 1-scctions and 2-scction are related to the canonical Liouville section of 
{jEge*^Y — ' E'g and the canonical symplectic section of a'^(T'^<^ E^)* — > E^ 
by pull-back under the discrete Legendre transformations, that is, 

{T¥-Ld,¥-LdreEa = ez,, {T¥+Ld,¥+LdreEa - e+_^, 

iT¥-Ld,¥-LdrnEa = {T¥+Ld,¥+LdrnEa = ^l,- 

Discrete regular Lagrangians. A discrete Lagrangian Ld ■ G — > M is said to be 
regular if the set of solutions of the discrete Euler-Lagrange equations is locally the 
graph of a diffeomorphism, that is, there exists locally a unique discrete Lagrangian 
evolution operator : G — > G for Ld- In such a case, is called the discrete 
Euler-Lagrange evolution operator. In (see Theorem 4.13 in [3S]), we obtained 
some necessary and sufficient conditions for a discrete Lagrangian on a Lie groupoid 
G to be regular that we summarize as follows: 

Ld is regular •<==^ The Legendre transformation ¥^ Ld is a local diffeomorphism 
■<==^ The Legendre transformation ¥~ Ld is a local diffeomorphism 
<;=^ The Poincare-Cartan 2-section il^^ is symplectic 

G 

on the Lie algebroid Eq x Eq = V (3 ©g Va — > G. 

Locally, we deduce that Ld is regular if and only if for every local basis {Xi] of 
Sec(T) the local matrix {X i{X jLd)) is regular. 
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Discrete Hamiltonian evolution operator. If : G — > R is a regular La- 
grangian, then pushing forward to Eq with the discrete Legendre transformations, 
we obtain the discrete Hamiltonian evolution operator, : Eq — > E^ 
which is given by 

a, -F±L<ioa,o(F±Ld)-^ (6.10) 
Defining the prolongation T : T^^ E^ — > T^^E^ of by 

we deduce that (see 

Noether's theorem. In Discrete Mechanics is also possible to relate invariance 
of the discrete Lagrangian under some transformation group with the existence of 
constants of the motion. In fact, we will say that a section X of Eq is a Noether's 
symmetry of the Lagrangian Ld if there exists a function / G C°° (M) such that 

In the particular case when dLd{X'^^''^^) = ~XLd + XLd — 0, we will say that X 
is an infinitesimal symmetry of the discrete Lagrangian Ld- 

li Ld : G — > R is a regular discrete Lagrangian, by a constant of the mo- 
tion we mean a function F invariant under the discrete Euler-Lagrangc evolution 
operator ^l^, that is. F o = F. Then, we have the following result. 

Theorem 6.2 (Discrete Noether's theorem). [33] If X is a Noether symmetry of 
a discrete Lagrangian Ld, then the function F = QJ^^{X^^'^^) — a* f is a constant 
of the motion for the discrete dynamics defined by Ld- 

7. Classical Field Theory on Lie algebroids 

In this section, we study Classical Field Theories on Lie algebroids. We consider 
a fiber bundle : M — > N, a Lie algebroid structure on a vector bundle t|[ : E — > 
M and a surjective morphism of Lie algebroids tt: E — > TN over v. The physical 
interpretation of the above data is as follows: we will consider a field theory in 
which the fields arc the sections of the bundle v and the partial derivatives of the 
fields are parameterized by linear sections of tt. 

We will find the equations for the extremals of a variational problem which 
roughly speaking is the following: given a Lagrangian function L defined on the 
set of sections of tt, and a volume form lo on the manifold N , we look for those 
morphisms of Lie algebroids which are critical points of the action functional 

5($) = / L($)w. 

J N 

This is a constrained variational problem, because we are restricting the fields to 
be morphisms of Lie algebroids, which is a condition on the derivatives of $. 

7.1. Jets. We consider two vector bundles : E — > M and t^: F — > N and 
a surjective vector bundle map tt: E — > F over the map v. M — > N. Moreover, 
we will assume that v. M — > is a smooth fiber bundle. We will denote by 
K — > M the kernel of the map tt, which is a vector bundle over M . Given a point 
m G M, if we denote n = i>{m), we have the following exact sequence — > Km — > 
E,n — > Fn — > 0, and we can consider the set 3,„7r of splittings 4> of such sequence. 
More concretely, we define the following sets iL,„7r = {w: Fn — > Em | w is linear }, 
3mT^ = {</>£ i^jriTT \ 77 o (j) = idp^ } and V„i7r = {ip £ iljiiTT I TT o = }. Therefore 
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UrnT^ is a vector space, V^tt is a vector subspace of HmT^ and 3mT^ is an afSne 
subspace of -CmTr modeled on the vector space V^Tr. By taking the union, = 
UmeA/'C'mTr, Stt = U„iej\/3„i7r and Vtt = UmeM^mTr, we get the vector bundle 
TTio : -Ctt — > M and the affinc subbundle ttiq : ^tt — > M modeled on the vector 
bundle TTio : Vtt — > M. We will also consider the projection tti : ^tt — > -/V defined 
by composition tti = ly o ttiq. An element of tJmTr will be simply called a jet at the 
point m G A/ and accordingly the bundle ^tt is said to be the first jet bundle of tt. 

Notice that the standard case 122] is recovered when we have a bundle f : M — > 
N and one considers the standard Lie algebroids E = TM — > M and F = TN — > 
N together with the differential of the projection n ^Tv. TM — > TN. With the 
standard notations, we have that J^iy = diTv). 

Local coordinates on ^tt are given as follows. We consider local coordinates 
(x') on A'^ and (x^u^) on M adapted to the projection v. We also consider local 
basis of sections {ca] of F and {eQ,eQ,} of E adapted to the projection tt, that is 
TT o gf, = Cq o 1/ and TT o = 0. In this way {ec} is a base of sections of K . An 
element w in £jm7r is of the form w = {w\eb + w^Ca) ® e"', and it is in 3m tt if and 
only if w'^ = S^, i.e., an element (j) in Jtt is of the form = (ca + d) e°. If we 

set = (j)", we have adapted local coordinates {x^,u^,y") on ^tt. Similarly, 

an element ijj € V^Tr is of the form tp = ipa^a ® e°. If we set Uai'i') = "0" i have 
adapted local coordinates (x', u"*, y") on Vtt. As usual, we use the same name for 
the coordinates in an affinc bundle and in the associated vector bundle. 

An clement z G defines an affinc function z on tJm'"' by contraction z{(j)) = 
(z, <j)) where (•, •) is the pairing (z, 0) = Tr(zo(/)) = Tr(0oz). Therefore, a section 9 of 
iL*7r defines a fibcrwisc affinc function 9 on ^Itt, 0(0) = (^irio(i/>)i '?^) — Tr(07rio(0) ° 4')- 
In local coordinates, a section of L*tt is of the form 9 = (6*^ (a;)e^ + 6'J^(x)e") Bq, 
and the affinc function defined by 6* is 61 = 9'^{x) + 9^{x)y". 

Anchor. Consider now anchored structures on the bundles E and F, that is. we 
have two vector bundle maps pf'- F — > TN and pe-E — > TM over the identity 
in A^ and M respectively. We will assume that the map tt is admissible, that is 
Pp OTT = Tv o pe- Therefore we have 

Q p^(ea) = Pajr- + Pa 



PHea) = pL^ and % 

PE[ea) = Pa- 



with pI = pI{x), = Pa{x,u) and = 

The anchor allows us to define the concept of total derivative of a function with 
respect to a section. Given a section a € Sec(i^), the total derivative of a function 
/ S C°°(M) with respect to a is the function df®a, i.e., the affine function 
associated to c?/ cr £ Sec(£;*7r). In particular, the total derivative with respect to 
an element Ca of the local basis of sections of F, will be denoted by f\a- In this 
way, if tr = a"'ea then df ® a = f\aO'°', where the coordinate expression of f^^ is 

Notice that, for a function / in the base A^, we have that f^a = P]i§§^ are just the 
components of df in the basis {eo}. 

Bracket. Finally, let us assume that we have Lie algebroid structures on rf^ : F — > 
N and on : E — > A/, and that the projection tt is a morphism of Lie algebroids. 
This condition implies the vanishing of some structure functions. 
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We have the following expressions for the brackets of elements in the basis of 
sections 

r Iea,ef,l = C;^fce^ + C^fcCc 
Iea,ef,] = C°f,ec and < l^a, ef3l= Cj^^e^ 

[ lea,ef3j = Clfje^, 

where C^^ = C^c(^) is a basic function. 

The structure functions can be conveniently combined as terms of some afhne 
function as follows 

In particular, such functions will appear in the Euler-Lagrange equations of the 
variational problem. 

7.2. Morphisms and admissible maps. By a section of tt we mean a vector 
bundle map <i> such that tt o $ = idp, (i.e., we consider only linear sections of tt 
(see also ^ni))- It follows that the base map N — > M is a section of i^, i.e., 
V o (f) = iAn- The set of sections of tt will be denoted by Sec(7r). The set of those 
sections of tt which are a morphism of Lie algebroids will be denoted by M(7r). We 
will find local conditions for $ 6 Sec(7r) to be an admissible map between anchored 
vector bundles and also local conditions for <I> to be a morphism of Lie algebroids. 

Taking adapted local coordinates (a;*,u'^) on Af, the map (f) has the expression 
4>{x^) = {x^ tU^{x)). If we moreover take an adapted basis {ca, Cq} of local sections 
of E, then the expression of $ is given by ^{ca) = Ca + y2{x)ea, so that the map $ 
is determined by the functions (u^{x), y"{x)) locally defined on N. The action on 
the dual basis is (E>*e° — e'^, and <I>*e" = y^(x)e°, and for the coordinate functions 
= and ^*u^ = ■u^(x). 

The admissibility condition reads ^*{df) ~ d{^* f) for every function / G 
C°°(M). Taking / ~ a;* we get an identity, while taking f = wc get the 
condition ^ 

In addition to the admissibility condition, the morphism condition reads ^*dd = 
d{^*0) for every section 9 oi E*. For = e° we get an identity, while for 6 = 
we find the 

i ^Vc i (^Vb _ „.afia , fia ,,0„,7 i ,,7 _ r'" A- — D 

dx'- dx' ^PlVb Vc + ^b-iVc '■^cjUb + '^bc — 

7.3. Variational Calculus. In what follows in this paper we consider the case 
where the Lie algebroid F is the tangent bundle F = TN with pp = idrw and [•, •] 
the usual Lie bracket of vector fields on N. The Lie algebroid E remains a general 
Lie algebroid. Moreover, for local expressions on F, the local basis of sections of 
F which we will consider is a basis of coordinate vector fields = so that 
Pl = SI and C," = 0. 

Variational problem. Given a Lagrangian function L G C°°{3tt) and a volume 
form u! G /\^ (TN), where r = dim(iV), we consider the following variational prob- 
lem: find the critical points of the action functional S{^) = /^L($)a; defined on 
the set of sections of tt which are moreover morphisms of Lie algebroids, that is, 
defined on the set M(7r). Here by i($) we mean the function n i-^ L($„), where 
$n G St^ is the restriction of $ the fiber Fn = TnN . 

It is important to notice that the above variational problem is a constrained 
problem, not only because the condition tt o $ = idj?, which can be easily solved. 
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but because of the condition of $ being a morphism of Lie algebroids, which is a 
condition on the derivatives of <i>. Taking coordinates on N such that the volume 
form is = dx^ A • • • A dx^\ the problem is to find the critical points of 

L{x\u'^{x),y^{x))dx^ A • • • Adx'', 



IN 

subject to the constraints 

— i n^ii" nnH ^^'^ i r"" i,'^ ii'^ 4- r° ifiiH J- r*" — n 

~ Pa + PaVa ^Ud - + C-b^l/c ~ L'c7% + '-^PiVb Vc + '-'be ~ U. 

The first method one can try to solve the problem is to use Lagrange multipliers. 
Nevertheless, one has no warranties that all solutions to this problem are normal 
(i.e., not strictly abnormal). In fact, in simple cases such as the problem of a heavy 
top PHIj one can easily see that there will be strictly abnormal solutions. Therefore 
we take another approach, which consists of finding explicitly finite variations of a 
solution, that is, defining a curve in M(7r) starting at the given solution. 

Variations and infinitesimal variations. In order to find admissible variations, 
we consider sections of E and the associated fiow. With the help of this flow we 
can transform morphisms of Lie algebroids into morphisms of Lie algebroids. 
Flov^r defined by a section. We recall that every section of a Lie algebroid has 
an associated local flow composed of morphisms of Lie algebroids 021 IM] ■ More 
explicitly, given a section cr of a Lie algebroid E, there exists a local fiow $<, : E — > 
E such that 

as s=o 

for any section 9 oi /\ E. Moreover, for every fixed s, the map 4>s is a morphism of 
Lie algebroids, and the base map (f>s : M — > M, the (ordinary) fiow of the vector 
field p{a) e X(M). 

Complete lift of a section. In this section we will define the lift of a projectable 
section of i? to a vector field on ^tt, in a similar way to the definition of the first jet 
prolongations of a projectable vector field in the standard theory of jet bundles |52| . 

We consider a section ct of a Lie algebroid E projectable over a section t] of F. 
We denote by the fiow on E associated to a and by <&s the flow on F associated 
to rj. We recall that, for every fixed s, the maps "ifs and $s are morphisms of Lie 
algebroids. Moreover, the base maps tps and (j)s, are but the fiows of the vector 
fields pE{<y) and pf{ti), respectively. 

The projectability of the section implies the projectability of the fiow. It follows 
that (locally, in the domain of the fiows) we have defined a map iL'I's : Lit — > 
by means of L'^ s{w) = 5',, o w o By restriction of L"^ s to ^Itt wc get a map 

which is a local fiow in ^tt. We will denote by X^^' the vector field on ^tt 
generating the fiow The vector field X^^^ will be called the complete lift to 

tJTT of the section a. Since tJ^'^ projects to the fiow ips it follows that the vector 
field X'-J-^ projects to the vector field pe{<^) in M. 

Locally, a section a — a'^Ca + cr"eQ, is projectable if a°- = (T°(a;') depends only 
on a;\ Its complete lift has the local expression 

where cr^^ = (t[^ + Z^^cr^ + Z^^^a^ - (^a^^ + cr^C^,,^ . In particular, if a projects to 
the zero section, i.e., a" ~ 0, wc have 
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Euler-Lagrange equations. Let $ S M(7r) be a critical point of S. In order 
to find admissible variations wc consider a 7r-vertical section a of E. Its flow 
■ E — > E projects to the identity in F = TN. Moreover we will require a to 
have compact support. Since for every fixed s, the map ^'s is a morphism of Lie 
algcbroids, it follows that the map = 4'^ o $ is a section of tt and a morphism 
of Lie algebroids, that is. s ^ is a curve in M(7r). Using this kind of variations 
we have the following result. 

Theorem 7.1. |42j Select local coordinates such that the volume form is expressed 
as u! ~ dx^ A • ■ • A dx^ . A map ^ is a critical section of S if and only if the 
components y" of <I> satisfy the system of partial differential equations 

A An 



dx°- 



dVa dy^ 
dxb dx^ ^ ''^^^ 



d 

dx" 



dL 
dy^ 



dL 
dyl 



dL 



du^ 



Pt 



Proof. Recall that by L{^) we mean the function in N given by L{^){n) = L{^n) 
where is the restriction of $: F — > E to the fiber F„. The function L($s) is 

L($,)(n) = L(*, o $„) = L(a^',($„)) = (a*:L)($)(n), 

and therefore the variation of the action along the curve s <E>s is 







ds 



s=0 



N 



ds 



s=0 



(i:^a)L)($)w. 



N 



Taking into account the local expression of for a 7r-vertical a, we have that 



dL 
du^ 
^^dL_ 

"du^ 



da^ 

dx'' 



dL 
dyl 





dL 


dyS 


d 


/ dL 


dx'' 


\dy2 



d 

dx'^ 



,dL_ 



Let us denote by 5L the expression with components 



Pt 



dL 



dyl "^du^' 

^-^^uJa with LUa ~ i _a_uj . Then we 



and by the (r — l)-form (along tti) J„ = 
have that 

= ^S{^s) = - / {5Ln + / d{J, o $). 
ds s=o Jf^ Jjv 

Since a has compact support the second term vanishes by Stokes theorem. More- 
over, since the section a is arbitrary, by the fundamental theorem of the Calculus 
of Variations, we get (5i = 0, which arc the Euler-Lagrange equations. Notice 
that the first two equations in the above statement are nothing but the morphism 
conditions. □ 



7.4. Examples. 



Standard case. In the standard case, we consider a bundle v. M — > iV, the stan- 
dard Lie algcbroids F = TN and E = TM and the tangent map n ^Tv: TM — > 
TN. Then we have that ^tt ~ J^v. When wc choose coordinate basis of vector 
fields (i.e., of sections of TN and TM) we recover the equations for the standard 
first-order field theory. Moreover, if we consider a different basis, what we get are 
the equations for a first-order field theory written in pseudo-coordinates [SJ 021 ■ 
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Time- dependent Mechanics. In [501 1511 we developed a theory of Lagrangian Me- 
chanics for time dependent systems defined on Lie algebroids, where the base man- 
ifold is fibered over the real line R. Since time-dependent Mechanics is nothing but 
a 1-dimensional field theory, our results must reduce to that. 

The morphism condition is just the admissibility condition so that, if we write 
x'^ = t and Uq = y", the Euler-Lagrange equations are 

dw^ A An, d f dL\ dL fs^ dL a 

in full agreement with |44| . In particular, for an autonomous system on a Lie 
algebroid V — > Q, one considers E = TR x V — > M x M with tt the projection onto 
the first factor TR. Our results provide yet another indication of the variational 
character of autonomous mechanical systems on Lie algebroids. 

Topological field theories. Given a closed r-form f2 on a Lie algebroid V — > Q, 
we can define a topological field theory as follows. For an r-dimensional manifold 
N we consider F = TN — > N, E = TN xV — > N x Q and tt the projection 
onto the first factor TN. The Lagrangian of the theory is L{^) = Then it is 

easy to see that the Euler-Lagrange equations reduce to the morphism condition. 
In this way, one can cope with systems such as Poisson cr-models or Chern-Simons 
theories 

Systems with symmetry. The case of a system with symmetry is very important 
in Physics. We consider a principal bundle v : P — > M with structure group G 
and we set = M, F = TN and E = TP/G (the Atiyah algebroid of P), with 
7r([u]) = Tv{v). Sections of tt are just principal connections on P and a section is 
a morphism if and only if it is a flat connection. The kernel K is just the adjoint 
bundle {Pxq)/G — * M . By an adequate choice of a local basis of sections of F, K 
and E one easily find the covariant Eulcr-Poincare equations [3 IHj . The covariant 
Lagrangc-Poincarc equations [Hj can also be recovered within this formalism. 

8. Future work 

We have illustrated the generality of the theory of Lie algebroids and groupoids 
in a wide range of situations, from nonholonomic Lagrangian and Hamiltonian 
systems and mechanical control systems to Discrete Mechanics and extensions to 
Field Theory. Current and future directions of research include the following: 

Hamilton- J acobi equation for a Hamiltonian system on a Lie algebroid: It wo- 
uld be interesting to continue with the study started in |26j of Hamilton- 
Jacobi theory for Hamiltonian systems on Lie algebroids. In particular, it 
would be interesting to introduce a suitable definition of a local (global) 
complete integral of the Hamilton-Jacobi equation. The idea would be 
that the knowledge of an integral of the equation would allow the "direct 
determination" of some integral curves of the corresponding Hamiltonian 
vector field. 

Geometric formalism for Vakonomic Mechanics on Lie algebroids: An inter- 
esting topic to study is the case of constrained variational problems on Lie 
algebroids. In this case, to derive the equations of motion for a Lagrangian 
system subject to nonholonomic constraints, one invokes a variational prin- 
ciple, rather than the Lagrange-D'Alembert's principle (cf. Section^)). The 
differential equations obtained, called vakonomic equations, are in general 
different. From an optimal control perspective, it seems interesting to 
generalize the formalism developed in |11|. 
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Mechanical control systems on Lie algebroids: Topics of interest related to 
mechanical control systems on Lie algebroids include the investigation of 
controllability tests along relative equilibria and the study of systems that 
include gyroscopic and dissipative forces. 

Discrete Mechanics on Lie groupoids: We are currently studying the construc- 
tion of geometric integrators for mechanical systems on Lie algebroids. 
We have introduced the exact discrete Lagrangian in the Lie groupoid 
formalism, and arc discussing different types of discretizations of continu- 
ous Lagrangians and their numerical implementation. We plan to explore 
natural extensions to forced systems and to systems with holonomic and 
nonholonomic constraints as in ^1 . 

Classical Field Theory and Lie algebroids: In |16[ I44| the authors have intro- 
duced the notion of a Lie affgebroid structure |17l 1411 HI) (see also |21|V 
They have developed a Lagrangian (and Hamiltonian) formalism on Lie af- 
fgebroidS; which generalizes some classical formalisms for time-dependent 
Mechanics and, in addition, may be applied to other situations. Since 
time-dependent Mechanics is a 1-dimensional field theory, it would be in- 
teresting to define the notion of a "Lie multialgebroid" , as a generalization 
of the notion of a Lie affgebroid. This mathematical object should en- 
code the geometric structure necessary to develop field theories. The first 
example of a Lie multialgebroidj should be ^Itt. The notion of a Lie multi- 
algebroid will potentially allow to study other aspects of the theory, such 
as Tulczyjew's triples associated with a Lie multialgebroid and Hamilton- 
Jacobi equation for classical field theories on Lie multialgebroids. 
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